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A  series  of  tests  of  incremental  stress  waves  of 

small  amplitude  were  performed  on  annealed  aluminum  1100-F 

rods  and  copper  rods  under  axial  quasi  static  preloading 

beyond  the  elastic  limit.   Particle  velocity  records  of  the 

incremental  waves  were  obtained  at  four  gage  stations  along 

the  rod.   The  strain  rates  of  the  quasi  static  preloading 

-6       -4     -1 
were  in  the  range  of  10   to  10    sec   ,  while  the  strain 

rates  of  the  order  of  1.0  sec~   were  achieved  in  the  incre- 
mental plastic  waves  with  maximum  particle  velocity  around 
50  cm/sec  and  duration  of  around  500  ysec. 

Test  results  indicated  that  except  near  the  plateau 
in  the  records  for  aluminum  any  given  level  of  incremental 
particle  velocity  propagates  along  the  rod  with  a  constant 
speed,  not  affected  by  the  strain  rate  within  the  small 
range  of  the  incremental  wave  strain  rates  encountered. 


IX 


However,  the  obtained  propagation  speed  function  for  each 
incretnental  wave  was  found  dependent  on  both  prestrain  and 
prestrain-rate  values.   A  precondition  parameter  composed 
of  prestrain  and  prestrain  rate  was  proposed  to  give  a 
quantitative  description  of  this  phenomenon.   This  discovery 
revealed  another  physical  property  of  metal  under  incremental 
loading  in  addition  to  the  well-known  phenomenon  of  the 
leading  edge  propagating  at  the  elastic  bar-wave  speed. 

Both  strain-rate- independent  and  strain-rate- 
dependent  theories  were  used  to  interpret  the  experimental 
results.   By  applying  the  simple-wave  solution  of  a  modified 
rate- independent  theory,  incremental  stress-strain  curves 
were  calculated  from  the  wave-propagation-speed  function  for 
both  cases  of  aluminum  and  copper.   Quasilinear  and  semi- 
linear  rate-type  constitutive  models  with  linearized  terms 
were  proposed.   Numerical  computer  solutions  were  obtained 
for  incremental  waves  at  various  prestrain  levels  of  both 
aluminum  and  copper  specimens.   Results  show  that  the 
proposed  quasilinear  constitutive  equation  is  adequate  for 
copper  and  the  semilinear  constitutive  equation  is  qualita- 
tively  adequate  for  aluminum.   Results  for  copper  showed 
that  matching  the  particle  velocity  records  by  adjusting 
the  parameters  in  the  quasilinear  constitutive  equation 
does  not  guarantee  a  reasonable  solution  in  incremental  stress 
and  strain,  when  it  was  compared  with  the  incremental 


stress-strain  curve  of  rate-independent  theory.   Therefore 
a  strong  inequality  was  proposed  for  determining  the  para- 
meter in  the  instantaneous  response  term  of  the  quasilinear 
equation.   This  proved  to  be  an  efficient  aid  in  the  para- 
meter identification  since  only  one  parameter  was  left, 
which  can  be  determined  from  matching  the  velocity  records. 
Results  show  that  the  parameters  in  the  rate-type  constit- 
utive equation  in  both  materials  depend  on  the  proposed 
precondition  parameter. 
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CHAPTER   1 
INTRODUCTION 

1  .  1   Background 

Dynamic  plasticity  is  the  branch  of  the  mechanics  of 
solids  in  which  rapid  deformation  of  solid  bodies  beyond  the 
elastic  limit  is  studied  without  neglecting  inertia  effects. 
The  transient  spreading  out  of  the  deformation  from  the 
loaded  part  of  a  body  subjected  to  impact  loading  is  refer- 
red to  as  plastic  wave  propagation.   When  inertia  effects 
are   included,  the  analysis  of  the  deformation  of  a  rapidly 
loaded  body  requires  the  consideration  of  plastic  wave  pro- 
pagation except  for  very    small  specimens  in  which  the  as- 
sumption of  uniform  stress  and  strain  may  be  a  satisfactory 
approximation  for  practical  purposes.   The  analysis  of  plas- 
tic wave  propagation  is  relatively  new,  with  a  history  of  about 
35  years.   There  was  some  study  of  the  subject  before  1940; 
the  first  important  paper  was  by  Donnell  (1930). 

Two  basic  assumptions  about  the  constitutive  law  of 
the  material  behavior  under  dynamic  plastic  loading  were 
proposed  in  the  1940's.   A  stra i n-ra te- i ndependent  theory 
was  proposed  independently  by  Taylor  (1940)  in  England, 
von  Karman  (1942)  in  America,  and  Rakhmatulin  in  the  Soviet 
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Union.  They  analyzed  the  problem  of  longitudinal  impact  of 
thin  rods  and  assumed  that  radial  inertia  effects  were  neg- 
ligible. A  si ngl e- val ued  stress-strain  relation,  independ- 
ent of  strain  rate  was  used  in  these  analyses.  The  partial 
differential  equation  of  motion  was  derived  and  a  simple 
wave  solution  for  continued  impact  on  a  semiinfinite  bar  was 
obtai  ned . 

Experiments  were  carried  out  shortly  after  the  devel- 
opment of  this  theory.   A  summary  paper  by  Duwez  and  Clark 
(1947)  presented  the  major  results  of  the  experimental  work 
in  the  preceding  few  years.   In  general,  the  prediction  of 
the  formation  of  a  strain  plateau  near  the  impact  end  in  a 
constant-velocity  impact  and  the  relation  between  the  impact 
velocity  and  the  maximum  plastic  strain  reached  were  fairly 
well  verified  by  the  experiments.   Some  systematic  discrep- 
ancies in  the  magnitude  of  the  residual  strain  and  in  the 
stress-time  history  at  the  fixed  end  opposite  the  impact  end 
were  observed.   In  the  tension  impact  tests  the  maximum  re- 
sidual strain  was  smaller  than  predicted  by  the  theory,  and 
the  observed  stress  time  variation  at  the  fixed  end  during 
impact  showed  the  stress  there  was   greater  than  the  theory 
predicted.   It  was  suggested  that  this  might  be  attributable 
to  strain-rate  effects  in  the  material  under  impact  loading. 

The  strain-rate-dependent  theory  was  proposed  by 
Sokolovsky  (1948)  and  by  Malvern  (1949)  independently  of 
each  other.   They  considered  the  plastic  flow  stress  as  a 


function  of  strain  rate  as  well  as  the  level  of  strain. 
Both  Sokolovsky  and  Malvern  proposed  semilinear  constitutive 
laws  and  used  them  to  analyze  examples  of  longitudinal  im- 
pact on  bars.   Mai vern ( 1 949  )  also  suggested  a  more  general 
quasi  linear  law  and  showed  how  it  could  be  used  to  study 
wave  propagation  by  the  method  of  characteristics.  Sokolovsky 
(1948)  assumed  perfect  plasticity,  while  Malvern  (1949)  con- 
sidered work  hardening.   The  constitutive  law  proposed  by 
Malvern  was  of  the  el  as ti c-vi sco-pl asti c  form  given  by 


Ee  =  o  +  g(a ,£ ) 


(1.1) 


where  g{o,e)    is  an  arbitrary  function  expressing  the  strain- 
rate  dependence.   An  idealized  form  of  the  constitutive  law 
(1.1)  was  used  in  which  g(a,e)  =  k[a  -  f(e)],  where  [o  -  f(e)] 
was  called  the  "overstress" ,  that  is,  the  excess  of  dynamic 
stress  over  the  static  yield  stress.   A  numerical  example 
was  given  using  the  idealized  linear  overstress  form,  and  it 
was  concluded  that  the  assumed  constitutive  law  could  ac- 
count for  discrepancies  observed  in  the  stress-time  varia- 
tion at  the  fi xed  end  of  an  impact  specimen.   However,  this 
type  of  flow  law  failed  to  predict  a  strain  plateau  near  the 
impact  end  of  constant-velocity  impact  specimens  in  the 
short  time  forwhich  the  calculations  were  made. 

The  question  concerning  the  possibility  of  the  exist- 
ence of  a  strain  plateau  in  the  solution  of  a  Malvern-type 
theory  was  considered  by  several  authors.   Rubin  (1954) 


obtained  a  solution  by  using  Laplace  transformation  for  the 
problem  of  the  small  increment  with  linearized  quasistatic 
curve.   An  asymptotic  plateau  was  proved  to  exist  from  the 
closed  form  solution.   However,  no  example  was  given;  how 
fast  the  asymptotic  plateau  was  approached  still  remained 
a  question.   Later  Bianchi  (1963),  Efron  (1964)  and  Wood  and 
Phillips  (1967)  showed  directly  the  asymptotic  plateau  in 
their  numerical  computer  solutions.   A  recent  theoretical 
work  by  Suliciu,  Malvern  and  Cristescu  (1972)  has  placed  re- 
strictions on  the  form  of  the  constitutive  equations  under 
which  a  plateau  can  occur  and  has  shown  that  with  the  usual 
form  of  the  Malvern  type  semilinear  models  only  an   asymp- 
totic plateau  is  possible. 

Concerning    the     apparent  increase  in  the  yield 
stress  at  high  velocity  impact,  Malvern  theory  predicts  such 
an  increase  in  the  value  of  initial  yield  stress  at  infinite- 
ly high  strain  rates  because  plastic  strain  does  not  occur 
instantaneously  but  takes  time  in  which  to  become  appreci- 
able.  At  high  finite  rates  the  dynamic  curves  rise  near  the 
elastic  line,  but  any  increase  in  the  yield  stress  would 
have  to  be  made  an  explicit  assumption.   Regarding  this 
point,  the  experiments  of  Campbell  (1954)  on  steel,  Tietz 
and  Dorn  (1949)  on  aluminum,  and  Taylor  (1957)  on  steel  all 
show  that  the  dynamic  yield  stress  increases  with  the  rate 
of  strain. 


Rubin  (1954)  also  showed  that  according  to  the  line- 
arized Malvern  theory  the  leading  edge  of  an  incremental 
plastic  wave  would  travel  at  the  elastic  wave  speed  instead 
of  at  the  plastic  wave  speed  predicted  by  the  rate-independent 
theories.   Experiments  by  Bell  (1951)  on  steel  bars  and 
Sternglass  and  Stuart  (1953)  with  copper  strips  involved 
the  propagation  of  incremental  impact  loads  superimposed 
upon  static  dead  loads  in  excess  of  the  elastic  limit.   The 
wave  fronts  were  found  to  propagate  at  the  elastic  wave 
speed.   Later  results  of  incremental  loading  waves  in  dy- 
namically prestressed  bars  were  obtained  by  Alter  and  Curtis 
(1956)  for  lead  and  by  Bell  and  Stein  (1962)  for  aluminum. 
These  results  also  showed  that  the  leading  incremental  wave 
edge  propagates  at  the  elastic  bar-wave  speed.   Similar  re- 
sults were  also  reported  by  Papirno  and  Gerard  (1961)  on  an- 
nealed aluminum  1100-H,  and  by  Bianchi  (1963)  on  copper 
strip. 

This  phenomenon  of  incremental  wave  was  widely  inter- 
preted as  evidence  of  strain-rate  effect  as  predicted  by  the 
Malvern  strain-rate-dependent  theory.   However,  Bell  and 
Stein  (1962)  disagreed  with  this  idea  and  stated  that  such 
an  interpretation  is  inconsistent  for  annealed  aluminum,  an- 
nealed copper,  and  pure  lead,  in  that,  large  amplitude  wave 
propagation  has  been  shown  to  be  in  very  close  agreement  with 
the  stra i n-ra te- i ndependent  finite  amplitude  theory  [see 
Bell  (1956),  (1960a),  (1961a, b,c),  (1968)].   They  explained 


the  phenomenon  by  a  trigger  mechanism  which  has  the  property 
of  rei naugurati ng  plastic  deformation  for  the  incremental 
wave  when  it  occurs  after  the  dynamic  prestressed  wave  has 
reached  a  maximum  value. 

Bell  and  his  coworkers  have  obtained  many  valuable 
experimental  results  in  finite  amplitude  stress  wave  propa- 
gation.  By  using  the  diffraction  grating  technique,  Bell 
(1960b)  has  shown  that  the  dynamic  overstress  is  determina- 
ble from  a  consideration  of  the  three  dimensional  develop- 
ment of  the  plastic  wave  in  the  first  diameter.   Also,  by 
means  of  diffraction  grating  measurements  of  strain  in  the 
free  flight  impact  of  identical  specimens  of  varying  dia- 
meter, the  geometric  character  of  the  wave  development  given 
in  Bell  (1961a)  was  shown  to  be  consistent.   Thus,  dynamic 
overstress  was  seen  to  depend  upon  the  diameter  and  not  upon 
the  distance  from  the  impact  face  as  would  be  expected  if 
the  phenomenon  were  associated  with  strain  rate  behavior. 
However,  Bell  maintains  that  the  dynamic  stress-strain  re- 
lation is  not  necessarily  the  same  as  the  quasi  static  rela- 
tion.  Similar  results  were  also  given  by  Kolsky  and  Douch 
(1962)  for  copper  and  aluminum. 

By  examining  various  forms  of  g(a,e)  in  equation  (1.1), 
Ripperger  and  Watson  (1968)  concluded  that  wave-propagation 
experiments  with  long  initially  unstressed  specimens  in  com- 
pression were  of  little  value  in  deciding  whether  dynamic 
material  behavior  obeys  the  rate- independent  or  the 


rate-dependent  theory.   Results  for  titanium  obtained  by 
Lawson  and  Nicholas  (1972)  in  torsion  also  support  this 
finding.   Recently  from  resultsof  an  experiment  which  could 
measure  simultaneously  tensile  and  compressive  wave  profiles 
in  the  specimen,  Khan  (1973)  demonstrated  the  strain-rate- 
independent  theory  is  also  applicable  for  propagation  of 
tensile  waves  in  annealed  aluminum  1100-F.   The  dynamic 
stress-strain  curve  was  found  to  be  the  same  in  tension  and 
compression. 

On  the  side  of  theoretical  development,  a  general 
quasi  linear  constitutive  equation,  which  includes  both 
the  rate- independent  theory  and  the  Malvern  semilinear  the- 
ory, as  special  cases,  was  proposed  by  Cristescu  (1963)  and 
Lubliner  (1964).   It  has  the  form 


£  -    ({)(a,£)a+  i|)(o,e) 


(1.2) 


with  the  implication  that  a  certain  part  of  the  plastic 
strain  is  also  developed  instantaneously.   Working  with 
Bell's  experimental  data,  Cristescu  (1972)  has  shown  that 
the  above  constitutive  equation,  with  special  forms  for  the 
functions  >^   and  i)    ,  can  fit  reasonably  well  data  of  the  pro- 
blem of  longitudinal  impact  between  two  thin  rods.   Lubliner 
and  Valathur  (1959)  studied  wave  propagation  in  constant 
velocity  impact  by  applying  a  similar  constitutive  equation 
with  4)  derived  from  a  Ramberg-Osgood  curve.   A  more  recent 
theoretical  work  by  Suliciu,  Malvern  and  Cristescu  (1974) 


has  proposed  some  restrictions  on  the  quasi! i near  material 
function.   It  was  shown  that  the  incremental  wave  behavior 
implies  that  at  least  one  continuous  transition  from  the 
quasi  linear  equation  to  a  semi  linear  equation  occurs  in  the 
neighborhood  of  the  relaxation  boundary,  under  the  assump- 
tion of  certain  continuity  properties  of  <p    and  ip    for  both 
increasing  and  decreasing  stress. 

The  split  Hopkinson  pressure  bar  introduced  by  Kolsky 
(1949)  has  been  widely  used  as  a  practical  method  in  obtain- 
ing the  dynamic  flow  properties  of  materials,  including 
strain  rate  effects  [see  Hauser,  Simmons  and  Dorn  (1961)  and 
Lindholm  (1964)].   This  technique  has  been  criticized  be- 
cause of  the  lateral  bulging  in  the  short  specimen  during 
compression  and  on  other  grounds  [see  De  Vault  (1965)  and 
Bell  (1966)]. 

Experiments  in  torsional  waves  were  sought  to  elim- 
inate any  complicating  effect  due  to  the  radial  inertia  that 
arises  in  longitudinal  impact.   Baker  and  Yew  (1966)  devel- 
oped a  torsional  split  Hopkinson  bar  and  compared  their  ex- 
perimental shear  strain  time  records  both  with  rate-independent 
theory,  using  the  quasi  static  stress-strain  curve,  and  with 
a  linear  overstress  rate-dependent  theory  as  used  by  Malvern, 
The  conclusion  was  that  the  rate-dependent  theory  gave  bet- 
ter agreement  qualitatively. 

Convery  and  Pugh  (1968)  studied  the  propagation  speed 
of  torsional  plastic  waves  in  metals  pretwisted  statically 


into  the  plastic  range.   A  new  method  was  developed  in  which 
a  tubular  test  specimen  together  with  a  concentric  bar  of  a 
brittle  material  was  twisted  slowly  such  that  at  some  time 
after  the  specimen  was  stressed  beyond  its  yield  the  brittle 
bar  broke  suddenly  and  transmitted  a  plastic  torsional 
stress  increment  along  the  specimen.   It  was  found  that  the 
velocity  of  propagation  both  in  copper  and  in  mild  steel  was 
the  same  as  the  elastic  shear  wave  velocity. 

Since  then,  experimental  results  of  this  kind  have 
appeared  frequently  in  the  literature.   Data  obtained  at 
higher  rates  of  strain  can  generally  be  classified  as  re- 
sults of  one  among  two  different  methods  of  testing.   In  the 
first,  specimens  are  exclusively  deformed  at  a  dynamic  rate 
during  the  test,  while  the  second  method  involves  preloading 
the  specimen  beyond  its  elastic  limit  at  a  chosen  rate  be- 
fore subjecting  it  to  an  incremental  disturbance.   In  both 
methods,  a  long  specimen  can  be  used  along  which  the  propa- 
gation of  the  stress  pulse  is  measured  and  analyzed;  alter- 
natively, with  a  short  specimen,  the  complete  history  of  the 
flow  stress,  as  related  to  the  strain-rate  history  involved, 
can  be  obtained. 

Reports  concerning  the  torsional  wave  propagation 
[see  Baker  and  Yew  (1966),  Yew  and  Richardson  (1969),  Lawson 
and  Nicholas  (1972),  Banerjee  and  Malvern  (1975)]  have  shown 
that  theoretical  predictions  of  experiments  involving  plas- 
tic wave  propagation  and  interaction  in  metals  are  not  very 
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sensitive  to  changes  in  the  rate-dependent  constitutive  equa- 
tion.  The  fact  that  most  of  the  plastic  deformation  in  the 
experiments  took  place  within  a  narrow  band  of  strain  rates 
explains  the  insensitivity  of  the  results  to  the  form  of  the 
constitutive  equation, as  long  as  the  stress  in  this  narrow 
band  of  strain  rates  does  not  vary  greatly  with  strain  rate. 

Therefore  most  experiments  seeking  to  determine   rate 
effects  were  performed  by  using  the  short  specimen  on  the 
torsional  split  Hopkinson  bar,  where  a  wider  range  of  strain 
rates  can  be  obtained.   These  experiments  can  be  divided 
into  two  categories:   (1)  approximately  constant  strain-rate 
test  over  a  broad  range  of  strain  rates  and  (2)  tests  with 
the  specimen  prestressed  at  a  quasi  static  or  dynamic  strain 
rate  and  then  subjected  to  an  incremental  disturbance.   Re- 
sults of  the  second  type  of  experiment  show  an  initial  elas- 
tic responce  to  the  incremental  loading  followed  by  strain 
hardening  with  increasing  deformation. 

An  example  of  such  an  incremental  stress-strain  curve 
is  given  here,  after  Frantz  and  Duffy  (1972).   Figure  1.1 
shows  the  constant  rate  curves  together  with  a  curve  ob- 
tained by  suddenly  changing  from  the  lower  rate  to  the  high- 
er rate.   Klepaczko  et  al .  (1974)  interpreted  this  kind  of 
stress-strain  behavior  by  including  a  strain-rate-history 
effect.   It  was  suggested  that  the  influence  of  strain  rate 
on  the  flow  stress  can  be  divided  into  two  parts.   The  ini- 
tial part  is  due  to  the  existing  work-hardened  structure  at 
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Figure  1.1   Torsional  Split  Hopkinson  Bar  Records  in 
Aluminum  [after  Frantz  and  Duffy  (1972)] 


the  strain  level,  while  the  second  is  associated  with  the 
formation  history  of  that  structure.   The  total  stress  dif- 
ference At  between  the  two  stress-strain  curves,  each  ob- 
tained at  a  constant  strain-rate,  can  be  divided  into  two 
parts.   The  first,  denoted  by  At  ,  is  probably  developed 
by  a  single  thermally  activated  dislocation  mechanism.   The 
second  part,  denoted  by  At,  ,  is  related  to  the  two  different 
structures  at  the  same  value  of  strain.   The  difference  be- 
tween these  two  structures  is  probably  due  to  a  dynamic  re- 
covery process  which  occurred  to  a  greater  degree  in  the 
preceding  deformation  at  the  lower  strain  rate. 
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Reports    of    similar  experimentson  copper  and  ti- 
tanium were  given  by  Eleiche  and  Campbell  (1974).   Smooth 
incremental  stress-strain  curves  were  obtained  for  the  ma- 
terial considered,  without  the  knee  reminiscent  of  an  upper 
yield  point  found  by  Frantz  and  Duffy.   Eleiche  and  Campbell 
used  a  fading  memory  to  explain  the  strain-rate  history  ef- 
fect.  Santosham  and  Ramsey  (1970)  made  an  incremental  wave 
test  on  annealed  copper  axially  prestressed  by  dead  loading. 
Incremental  stress-strain  curves  were  calculated  at  various 
prestrain  levels  by  applying  a  modified  strain-rate-independent 
theory. 

1 . 2   Outline  of  the  Present  Study 


To  find  the  nature  of  strain-rate  effects  in  mate- 
rials under  dynamic  loading  has  been  the  major  concern  in  the 
past  history  of  dynamic  plasticity.   As  the  results  from  the 
torsional  split  Hopkinson  bar  confirmed  the  old  results  in 
axial  impact  loading,  there  has  been  fairly  general  agree- 
ment that: 

(a)  Strain-rate  effect  plays  an  important  role  in  the 
incremental  waves. 

(b)  A  strictly  one-dimensional  plastic  stress  wave, 
with  most  of  the  deformation  occurring  within  a 
narrow  band  of  strain  rates  as  in  most  of  the 
cases  considered,  can  be  predicted  well  either  by 
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the  strain-rate-independent  theory  with  a  single 
dynamic  stress-strain  curve  or  by  a  strain-rate- 
dependent  theory  with  a  suitable  choice  of  con- 
stitutive equations. 
The  study  reported  in  this  thesis  contains: 
i)   Experimental  results  of  incremental  wave  tests  on 
annealed  copper  and  aluminum  rods  under  quasistat- 
ic  preloading.   Incremental  waves  were  recorded  by 
measurements  of  particle  velocity  at  four  stations 
along  the  rod.   It  was  found  that  an  exponential 
function  of  two  parameters  can  fit  well  the  com- 
puted wave-propagation  speed  function.   The  major 
parameter  of  the  wave-propagation  speed  function 
was  found  to  be  dependent  on  both  prestrain  and 
prestrain  rate.   This  discovery  provided  direct 
evidence  that  the  incremental  waves  are  both  pre- 
strain and  prestrain-rate,  or  strain-rate-history, 
dependent . 
))   I  nterpretation  of  the  experimental  results  by  us- 
ing both  stra i n- rate-i ndependent  and  strain-rate- 
dependent  theories.   Incremental  stress-strain 
curves  were  obtained  from  the  fitted  function  of 
wave-propagation  speed  by  applying  a  modified 
rate- i ndependent  theory  similar  to  that  of Santosham 
and  Ramsey  (1970).   Numerical  computer  solutions 
were  performed  for  the  rate-dependent  theories  of 
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Malvern  and  Cristescu.   The  velocity  record  of  the 
first  gage  station  was  used  as  the  boundary  condi- 
tion to  predict  the  velocity  versus  time  at  the 
other  three  gage  stations.   A  new  criterion  was 
suggested  for  determining  adequate  parameter  val- 
ues to  get  a  good  solution. 
Chapter  2  describes  the  experiment.   Experimental  re- 
sults are  given  in  Chapter  3, and  an  exponential  function  is 
fitted  to  the  wave  speed  data.   Chapter  4  gives  theoretical 
analysis  and  numerical  solutions,  and  Chapter  5  presents 
discussion  and  conclusions. 


CHAPTER   2 
THE  EXPERIMENT 

2. 1    General  Description 

A  3-meter-high  steel  frame,  with  a  hydraulic  jack 
at  the  bottom  and  a  load  cell  at  the  top  was  built  for  the 
purpose  of  studying  incremental  stress  waves  in  a  prestressed 
rod.   A  series  of  tests  were  performed  on  several  aluminum 
and  copper  rods  of  0.95  cm  nominal  diameter. 

A  schematic  diagram  of  the  test  system  is  shown  in 
Figure  2.1,   Uniaxial  tensile  loading  was  achieved  by  using 
compressed  hydraulic  fluid  in  order  to  avoid  bringing  any 
kind  of  mechanical  vibration  to  the  system.   While  the 
specimen  rod  was  being  continuously  prestrained  into  plastic 
deformation  up  to  0.05  nominal  strain,  a  series  of  incremental 
waves  was  recorded  at  4  locations  along  the  rod  by  using 
electromagnetic  velocity  transducers.   The  transducers  are 
described  in  Section  2.3. 

The  incremental  particl e- vel oci ty  wave  generated  in 
an  elastic  material  was  found  to  be  approximately  a  square 
pulse  with  amplitude  around  100  cm/sec  and  rise  time  20  ysec. 
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Figure    2.1       Schematic    Diagram    of    Test    Set    Up 
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Unloading  occurs  when  the  first  compressive  wave  front  re- 
flected from  the  bottom  reaches  the  first  velocity  trans- 
ducer station.   The  loading  durations  of  the  incremental 
waves  were  about  300  ysec  and  400  ysec  for  aluminum  and 
copper  specimens,  respectively. 


2.2    Specimens 


Commercially  pure  copper  rods  and  aluminum  1100-F 
rods  were  chosen  as  specimen  materials.   Both  rods  were 
manufactured  by  cold  extruding  with  a  nominal  diameter 
0.95  cm  and  3.66  m  original  length.   The  specimens  of  2.05  m 
length  were  cut  down  from  them.   Four  holes  each  12.7  cm 
apart  were  drilled  at  the  velocity  transducer  locations  by 
using  a  high  speed  drilling  machine  with  Whitman  &  Barnes 
No.  72  drill  (0.0250  inch  diameter). 

Another  0 . 45-meter- 1 ong  rod  was  also  cut  down  from 
the  remaining  1.61  m  length  of  rod.   Both  long  and  short 
rods  were  then  annealed  simultaneously  and  furnace  cooled 
in  a  3-meter-long  electrical  furnace  which  was  built  in  the 
laboratory  [see  Figure  2.2].   Temperatures  and  times  used 
for  annealing  were  345+4°C,  0.5  hr  for  aluminum  and  370±4°C, 
2.0  hr  for  copper.   The  densities  of  the  two  kinds  of  mate- 
rial used  were  found  as 


Figure  2.2   Three-Meter-Long  Electric  Furnace  with 
a  High  Power  Proportional  Solid  State 
Temperature  Controller. 


al umi  num 
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2.67  ±  0.004  gm/cm" 


P.^r^n^v.  =  8.89  ±  0.003  gm/cm^ 
^copper  ^ 


2  .  3   Velocity  Transducer 


The  velocity  transducer  measures  the  voltage  produced 
in  a  transverse  conductor  attached  to  a  rod  of  nonmagnetic 
material  as  the  rod  moves  axially  in  a  magnetic  field.   See, 
for  example,  Ripperger  and  Yeakley  (1963)  and  Efron  and 
Malvern  (1969).   The  transducer  is  based  on  the  Faraday  Prin- 
ciple of  Electromagnet  ism,  which  can  be  expressed  as 


-  /  B  X  vdl 


(2.1) 


where  B  is  the  magnetic  field  intensity  vector,  v  is  the 
velocity  of  the  conductor  relative  to  the  magnetic  field,  and 
1  is  the  vector  representation  of  length  measured  along  the 
conductor.   The  units  in  the  SI  system  are 


e   :  volts 


t    :    weber/m   (or  10   Gauss) 
V  :  m/s 
T  :  m 


From  the  arrangement  shown  in  Figure  2.1,  it  can  be 
seen  that  the  longitudinal  displacement  of  each  velocity 
transducer  conductor  varies  as  much  as  4.5  to  6.5  cm  as  the 
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rod  is  carried  to  0.05  nominal  tensile  prestrain.   Therefore 
two  requirements  must  be  fulfilled  in  order  to  have  good 
velocity  transducers.   They  are    (1)   a  long  magnetic  field 
is  required  for  each  transducer,  and  (2)   the  velocity  trans- 
ducer must  be  able  to  move  along  with  the  rod. 

Four  pairs  of  permanent  horseshoe  magnets  were  ob- 
tained from  surplus.   The  original  area  of  each  magnetic 
pole  was  2.2  cm  x  6.0  cm  and  the  distance  between  north  and 
south  poles  was  2.8  cm.   The  area  of  each  magnetic  pole  was 
enlarged  to  4  cm  x  12  cm  by  adding  to  it  pieces  machined 
from  soft  magnetic  iron.   This  makes  the  magnetic  field  more 
uniform  and  thus  reduces  the  errors  from  possible  deviations 
in  position  measurement  of  the  velocity  transducer.   Figure 
2.3  shows  an  example  of  the  difference  in  the  magnetic  field 
intensity  before  and  after  the  enlargement  of  the  pole  area. 

Figure  2.4  shows  a  schematic  diagram  of  the  velocity 
transducer  which  was  designed  and  built  to  meet  the  above 
mentioned  second  requirement  after  several  other  kinds  of 
unsuccessful  attempts.   The  transducer  is  made  by  passing  a 
No.  25  gauge   enamel -coated  wire  through  a  hole  drilled 
along  the  diameter  of  the  rod  cross-section.   The  wires 
came  out  at  points  A'  and  D'  from  holes  on  the  brass  holder. 
It  is  by  the  inertia  effect  of  the  brass  holder  that  during 
the  passing  of  the  stress  wave  only  the  movement  of  section 
BC  is  detected.   The  total  displacement  of  section  BC  is 
0.4  mm  for  a  square  wave  of  100  cm/sec  amplitude  and  400  ysec 
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Figure    2.4      Incremental    Particle    Velocity    Transducer, 
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Figure  2.5   Trace  Records  for  Outputs  of  Strain  Gage  and 
Velocity  Transducer. 


24 


duration.   For  the  incremental  waves  in  this  experiment  the 
displacement  is  much  less  than  0.4  mm.   This  means  that  the 
angle  of  rotation  of  wire  AB  or  CD  with  respect  to  point  A 
or  D  is  less  than  0.08  radian  or  4.6  degrees.   Four  pairs 
of  such  brass  holders  were  made. 

The  validity  of  the  velocity  transducer  was  also 
checked  by  obtaining  a  signal  simultaneously  from  a  surface 
strain  gage  on  an  aluminum  rod  for  an  elastic  wave.  Resist- 
ance strain  gages  of  Micro-Measurement  Type  ED-DY-250  BF-500 
with  gage  factor  3.31  were  selected  for  the  test.   Because 
the  amplitude  of  the  elastic  wave  in  strain  is  only  about 
0.0001,  this  type  of  gage  with  a  resistance  of  500  ohms  was 
selected  in  order  to  give  better  resolution.   The  gages  were 
mounted  in  pairs,  diametrically  opposed  to  cancel  any  effects 
due  to  bending,  at  a  position  30  cm  below  the  lowest  velocity 
transducer  station.   The  gages  were  connected  in  series  and 
installed  as  one  arm  of  a  Wheatstone  bridge.   The  results 
are  shown  in  Figure  2.5.   The  upper  curve  is  the  strain  gage 
output,  the  lower  curves  are  velocity  transducer  outputs. 
It  is  obvious  that  the  shapes  of  these  two  kinds  of  signal 
agree  well.   The  magnitude  of  each  signal  was  also  checked. 
The  results  are 


Station  No.    Velocity  (cm/sec)    Strain 


71.5 
71.3 

71.8 
71.9 


142x10' 


146x10' 


c  =  velocity/strain 


+5 
5.03x10   cm/sec 


4.92x10"^^  cm/sec 
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The  average  of  these  two  values  of  c   is  about  VL    below  the 
^  0 

value  of  elastic  bar  wave  speed  obtained  in  Section  3.3. 


2 . 4   Velocity  Transducer  Calibration 


An  aluminum  frame  of  length  0.5  m  which  can  span 
all    the    four  magnets  was  made  and  mounted  as  shown  in 
Figure  2.6.   Since  the  positions  of  velocity  transducers 
varied  as  the  rod  was  being  prestrained  and  the  magnetic 
field  was  not  uniform,  the  calibration  of  the  velocity  trans- 
ducer was  required  to  give  a  relation  between  the  effective 
magnetic  field  intensity  at  the  conductor  wire  and  the  con- 
ductor position  expressed  in  a  coordinate  system.   The  alu- 
minum frame  functioned  essentially  as  a  reference  coordinate 
frame,  on  which  the  magnetic  field  was  mapped  and  with 
which, later  in  each  series  of  incremental  wave  tests, the 
initial  and  final  positions  of  the  transducer  conductor  were 
measured . 

An   0.  S.  Walker   Company  type  MG-1  Gaussmeter  with 
Hall-effect  flat  probe  was  used  to  map  the  magnetic  field. 
The  coordinate  system  was  chosen  so  that  the  origin  was 
close  to  the  center  of  rod.   Field  intensity,  denoted  by  B  , 
was  measured  at  fifteen  points  on  an  x-y  plane  for  every 
0.5  cm  increment  of  z-coordinate.   In  practice,  (x,y)  was 
fixed  first,  and  the  probe  moved  along  z-axis,  stopping  for 
a  few  seconds  at  qnq.t)i    0.5  cm  interval.   The  output  from  the 


26 


Figure  2.6   Magnetic  Field  Mapping 
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Figure  2.7   Interpolation  Grids  and  Transducer  Positions 
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Gaussmeter  was  recorded  by  a  Honeywell  El ectroni k- 1 94  strip- 
chart  recorder  with  10  inch  pen  movement.   Each  stop  at  a 
0.5  cm  interval  showed  on  the  records  as  a  flat  step.   All 
the  data  were  read  and  punched  as  input  data  cards  for  the 
Velocity  Data-Processing  Computer  Program   [see  Section  3.3]. 

Because  the  transducer  wires  were  enclosed  by  the 
brass  holders,  so  that  their  positions  could  not  be  observed 
during  the  test,  five  reference  marks  were  made  on  the  rod 
at  known  positions  relative  to  the  wire  holes  before  the 
brass  holder  was  mounted  [see  Figure  2.7(c)].   During  the 
incremental  wave  test,  the  (x,y,z)  coordinate  values  of 
these  five  marked  points  were  measured  at  two  or  three  pre- 
strain  levels,  usually  the  beginning  and  end  of  each  dynamic 
incremental  loading.   Thus,  for  any  value  of  prestrain,  the 
(x,y,z)  coordinate  values  of  the  transducers  could  be  ob- 
tained by  interpolation.   Once  a  transducer ' s  coordinates 
are  determined,  denoted  by  (x^,  y^ ,  z,),  the  effective  mag- 
netic field  intensity  for  this  velocity  transducer  can  be 
obtained  by  a  subroutine  contained  in  the  velocity  data 
processing  computer  program.   This  routine  essentially  con- 
tains the  following  procedures: 


(a)  Determine  the  values  of  B  at  the  15  points  on 
the  z  =  Z-,  plane  by  using  linear  interpolation 
from  the  two  adjacent  z-planes  where  values  of 
B^  were  known  [see  Figure  2.7(a)]. 
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b)   Assuming  that  the  conductor  wire  BC  is  perfectly 
parallel  with  y-axis,  evaluate  B   at  9  points  on 
line  BC.   Successive  points  are  0.125  cm  apart. 
A  subroutine  automatically  determines  the  element 
number  to  which  the  points  belong.   The  equation 
of  interpolation  for  a  generic  point  shown  in 
Figure  2.7b  is  simply: 


+  r-(B,-Bj  +  t.(B,-Bj  +  r.t.(B  -B,+B  -B  ,)     (2.2) 


'b  "a 
where 


'd  "a 


abed 


X   -  X 
_P a_ 


y  -  y 
-^p   -^a 


(2.2a,b) 


(c)   By  applying  the  trapezoidal  rule,  the  calibration 
factor  for  converting  voltage  to  velocity  is  ob- 
tained as: 


/  B  dl  = 
C   "" 


2(n-l)  ^^1  ^ 


;  +  2(B,  + 
n     ^  2 


'n-1 


].BC   (2.3) 


where  n  =  9  in  this  case. 

From  equation  (2.1),  it  is  obvious  that 


V  -    e^/g 


(2.4) 


2.5   Quasistatic  U  n  i  a  x  i  a  1  Tension  Test 


The  purpose  of  this  research  was  to  study  the  incre- 
mental wave  on  the  quasi stati cal ly  prestrained  rod.  Therefore 
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a  complete  and  accurate  record  of  the  quasi  static  stress- 
strain  relation  was  needed.   Since  prestrain  had  to  be  meas- 
ured at  some  other  location  apart  from  the  velocity  trans- 
ducers, the  uniformity  of  the  long  rod  had  to  be  checked  at 
first.   The  major  factor  that  influences  the  uniformity  is 
the  annealing  temperature  distribution  along  the  rod.   The 
long  furnace  was  made  by  joining  three  cores.   It  was  found 
that  the  temperature  was  lower  near  the  junctions.   The 
position  of  the  rod  in  the  furnace  was  arranged  as  shown  in 
Figure  2.8.   The  effect  of  annealing  temperature  on  the 
specimen  rod  was  determined  by  checking  the  stress-strain 
curve  for  the  5  sections  A,  B,  C,  D,  E,  cut  down  from  a  7-ft 
aluminum  rod  after  annealing.   The  results  show  that  there 
are  hardly  any  differences  in  stress-strain  relation  among 
A,  C,  and  E  or  between  B  and  D.   However,  the  stress-strain 
curve  for  B  and  D  consistently  lies  above  that  of  A,  C,  and 

E  with  a  difference  in  stress  value  varying  monotoni cal ly 

2  2 

from  1  MN/m   for  small  strain  to  0.3  MN/m   for  6^=  0.05. 

Therefore  the  prestrain  value  obtained  from  the  extensometer 
located  in  section  A  should  represent  very   well  the  local 
strain  at  the  velocity  transducers,  which  were  located  in 
section  C  . 

The  load  cell  on  the  main  test  frame  was  calibrated 
by  adding  step  by  step  ten  10-lb  weights  before  and  after 
each  set  of  incremental  wave  tests.   The  aluminum  specimens 
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need  a  1000  lb  force  to  load  up  to  0.05  nominal  strain, 
while  the  copper  specimens  need  a  1700  lb  force  to  do  the 
same  thing.   The  force  was  supplied  from  an  AIRCO  Nitro- 
gen tank  with  full  pressure   2200  psi.   It  was  found  that 
in  order  to  obtain  a  1000  lb  force  the  pressure  needed 
was  1150  psi.   The  force  obtained  is  proportional  to  the 
pressure . 

Prestrain  records  for  incremental  wave  test  were 
obtained  by  using  a  5-inch  gage  length  extensometer  and 
also  a  2- inch  gage  length  LVDT  (Linear  Variable  Differential 
Transformer)  which  was  a  part  of  the  recording  equipment  of 
the  Tinius  Olsen  U-Celtronic  testing  machine.   The  LVDT 
was  calibrated  by  comparison  with  the  readi ng 'from  the 
extensometer,  when  both  instruments  were  mounted  on  a  short 
rod  loaded  on  the  Tinius  Olsen  Machine.   The  LVDT  was  mounted 
inside  the  5  inch  range  of  extensometer ' s  gage  length.   The 
result  of  the  calibration  showed  that  an  initial  nonlinear 
relation  exists  until  strain  larger  than  0.006,  after  which 
the  output  from  the  LVDT  is  a  linear  function  of  strain. 

A  6-volt  lantern  battery  was  used  to  drive  the 
VJheatstone  bridge  circuit  of  the  load  cell.   Output  from  it 
was  first  amplified  by  a  Hewlett-Packard  HP-425A  DC  Micro- 
Volt-Ammeter,  then  fed  into  one  channel  of  a  Gould  Brush 
220  2-channel  stripchart  recorder.   Output  from  the  LVDT 
was  fed  into  the  other  channel.   Strains  from  the  reading 
of  the  extensometer  gage  dial  were  also  recorded  on  the 
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Brush  220  by  a  remote  push-button  event  marker.   The  resol  ution 
of  the  strain  reading  by  the  extensometer  is  0.0001.  During 
the  incremental  wave  test,  the  Brush  220  chart  speed  was  set 
at  25  mm/min.  and  sensitivity  was  set  at  10  mv/div  for  both 
channels,  which  had  50  divisions  full  range. 

2 . 6   Incremental  Wave  Test  and  Wave  Recording  Equipment 


From  the  results  of  some  preliminary  tests,  it  was 
found  that  the  characters  of  the  incremental  waves  are  sub- 
stantially different  among  three  kinds  of  loading  states, 
which  are:   (1)  constant  prestressed  state,  which  is  the 
case  that  the  force  exerted  on  the  jack  piston,  by  the  pres- 
sured hydraulic  fluid,  is  in  equilibrium  with  the  resistance 
force  in  the  rod,  (2)  constant  prestrained  state  which  is 
the  case  that  the  manual  control  valve  was  closed,  and 
(3)  continuous  quasi  static  loading  state.   Since  creep  and 
relaxation  phenomena  were  involved  in  the  constant  pre- 
stressed and  prestrained  states,  respectively,  and  their 
influences  on  material  structure  are  not  well  known,  no 
further  incremental  wave  tests  were  done  on  these  two  pre- 
loading states. 

The  incremental  wave  was  generated  by  striking  a 
steel  tubular  hammer  on  a  collar  attached  on  the  bottom  end 
of  the  hard  aluminum  transmitter  tube.   The  top  end  of  the 
transmitter  tube  was  attached  to  the  specimen  rod  by  a 
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4-piece  grip  collar.   This  collar  was  located  15  cm  below  the 
lowest  velocity  gage  station  [see  Figure  2.1].   The  steel 
hammer  has  a  mass  of  345  gm.  The  dropping  height  is  15.6  cm. 
This  resulted  in  a  free-fall  striking  velocity  of  174  cm/sec. 
This  striking  velocity  was  increased  by  connecting  the  steel 
hammer  and  bottom  collar  with  3  soft  tension  springs.   The 
spring  constant  was  found  to  be  7.74x10   dyne/cm.   This 
arrangement  increased  the  striking  velocity  to  351  cm/sec. 

Some  examples  of  the  velocity  pulses  generated  by 
this  kind  of  hammer  mechanism  on  elastic  rods  are  shown  in 
Figure  2.9.   It  was   found    that  the  transmitted  velocity 
in  the  rod  depends  on  the  tightness  condition  of  the  grip. 

For  the  reasons  cited  in  the  first  paragraph  of  this 
section,  a  compressed  nitrogen  tank  with  an  output  pressure 
regulator  is  a  necessary  arrangement.   By  slowly  turning  up 
the  output  pressure  on  the  regulator,  the  specimen  rod  can 
be  further  prestrained  continuously.   This  avoids  the  relaxa- 
tion of  the  specimen  in  the  usual  procedure  of  closing  the 
manual  valve  first  and  increasing  the  pressure  in  the  accu- 
mulator to  a  higher  value,  then  opening  the  valve  and  re- 
loading.  Figure  2.10  is  a  general  view  of  the  physical  set 
up  of  testing  equipment. 

Outputs  from  the  four  velocity  transducers  were  fed 
into  four  Tektronix  1A7A  high-gain  differential  amplifiers 
in  rack-mounted  Tektronix  127  pre-ampl i f i er  power  supplies. 
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No.    1 

-^    ^ .         0. 5    mv 

Station        _t 


0^^^*:         0.?niv 
Stati ons 


(a) 
Initial  State,  e^  =  0.00 
Average   v  =  22.5  +  0.5  cm/sec 


All        ' 
^l\.            0.5  mv 
Stations   _j 


(b) 
Residual  Strain  =  0.0483 
Average   v  =  47.6  +  0.4  cm/sec 


Figure  2.9   Incremental  Elastic  Waves  (From  Copper 
Specimen  No.  4 ) . 
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Figure  2.10   General  View  -  Load  Frames 
Electronics  and  Recording 
Equ  i  pment . 
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The  high  and  low  frequency  -3dB  points  were  set  at  1  MHZ  and 
DC,  respectively.   The  outputs  from  the  1A7A  units,  which 
were  approximately  0.25  volts  for  an  input  of  1  mv  when  1A7A 
input  setting  was  1  mv/cm,  were  then  connected  to  the  inputs 
of  two  Tektronix  lAl  dual-trace  plug-in  units  on  a  Tektronix 
Type  556  dual  beam  oscilloscope.   All  four  channels  of  the 
two  lAl  plug-in  units  were  set  at  0.2  v/cm.   The  whole  re- 
cording system  was  calibrated  by  inputting  a  square  wave 
signal  into  one  1A7A  differential  amplifier  at  a  time.   The 
magnitude  of  the  square  wave  signal  was  set  equal  to  the 
value  of  the  scale  variable  setting  on  the  1A7A.   The  gain 
adjustment  on  the  1A7A  was  then  adjusted  to  make  the  square 
wave  trace  appeared  on  the  osci  1 1  oscope  to  have  a  1  cm  height. 
This  was  checked  once  again  after  each  set  of  incremental 
wave  tests.   The  gain  in  voltage  of  the  system  is  50  dB 
(400  times)  for  input  setting  on  1A7A  at  0.5  mv/cm  and  60  dB 
(1000  times)  for  0.2  mv/cm. 

The  time  base  on  the  oscilloscope  was  calibrated  by 
inputting  a  sine  wave  from  a  Wavetek  Model  112  wave  generator 
simultaneously  into  each  channel  and  a  Monsanto  Model  1 OOA 
Counter-timer.   The  exact  frequency  of  the  sine  wave  was 
obtained  from  the  counter.   The  display  on  the  oscilloscope 
screen  was  recorded  by  using  a  Polaroid  camera.   The  results 
of  the  time  base  calibration  for  the  100  ysec/cm  and  50  Msec/fcm 
ranges  on  the  oscilloscope  are  obtained  as  shown: 


100  ysec/cm  range: 

upper  beam  :  102.51  ysec/cm 

lower  beam  :  100,00  ysec/cm 

50  y  sec/cm  range  : 

upper  beam  :  50.80  ysec/cm 

lower  beam  :  49.78  ysec/cm 


A  Polaroid  camera  was  used  to  make  a  permanent  re- 
cord of  each  single-sweep  trace  on  the  oscilloscope.   The 
shutter  speed  was  set  at  1  sec  and  a  remote  control  push- 
button was  used  to  open  the  shutter  just  before  dropping 
the  hammer.   The  oscilloscope  was  set  at  single  sweep  and 
external  trigger.   The  trigger  was  provided  by  a  piezo- 
electric crystal  element  mountedonthe  top  of  the  steel  hammer, 
This  arrangement  allows  for  the  delay  in  the  scope  sweep 
mechanism. 


CHAPTER   3 
EXPERIMENTAL  RESULTS  AND  DATA  REDUCTION 

3 . 1   Quasistatic  Stress-Strain  Curves  and  Their  Slopes 

The  quasistatic  stress -strai n  curves  were  obtained 
from  load-time  and  strain-time  histories  recorded  during  the 
prestraining  of  each  specimen  rod  in  the  incremental  wave 
tests.   The  results  for  7  annealed  aluminum  specimens  and 
4  annealed  copper  specimens  are  shown  in  Figures  3.1  and 
3.2,  respectively. 

It  can  be  seen  from  Figure  3.1  that  the  maxumum  vari- 
ation in  stresses  among  different  aluminum  specimens  is 

2 
about  6  MN/m   which  is  around  25%  of  the  yield  stress  value. 

Results  for  copper  specimens  shown  in  Figure  3.2  appear  to 
be  more  uniform  among  different  specimens. 

The  large  variation  in  stresses  among  the  aluminum 
specimens  is  believed  to  be  caused  by  the  unequal  amounts 
of  previous  work  hardening  of  the  specimens  during  the  ex- 
trusion process  of  manufacturing.   The  annealing  process 
used  in  this  experiment  can  only  recover  about  75%  of  pre- 
vious work  hardening.   No  attempt  was  made  to  average  these 
curves,  since  each  curve  should  describe  the  best  stress- 
strain  relation  for  the  individual  specimen- 
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The  slope  of  the  stress-strain  curve,  dS/de,  versus 
strain,  e,  was  obtained  numerically  for  each  individual 
stress-strain  curve.   The  stress  value  corresponding  to  each 
Ae  strain  increment  was  first  obtained  by  interval-by- 
interval  cubic  Lagrange  interpolation.   The  value  of  dS/de 
at  e  =  e-  was  then  obtained  by  calculating  i^-^-t    -   S._-,)/(2  Ae). 
The  results  for  aluminum  and  copper  specimens  are  shown  as 
the  lower  curves  in  Figures  3.1  and  3.2,  respectively. 
Ae  =  0.00025  was  used  in  these  calculations. 

It  will  be  shown  later  in  Section  3.3  that  the  wave- 
propagation  speed  as  a  function  of  incremental  particle  ve- 
locity, denoted  by  c(v),  is  affected  by  both  prestrain  and 
prestrain  rate.   The  prestrain  rate  was  determined  by  two 
approaches  and  the  average  value  of  them  was  adopted.   The 
first  way  was  to  measure  the  slope  directly  from  e-t  record. 
In  practice  the  angle  of  slope  is  in  the  range  of  0.5  degree 
to  10  degrees.   The  error  in  reading  the  angles  is  large  in 
those  angles  less  than  5  degrees.   The  second  way  was  to 
measure  dS/dt  from  load-time  record  first,  then  to  obtain 
de/dt  from  de/dt  =  (dS/dt)/(dS/de).   Some  typical  examples 
of  the  values  of  de/dt  are  tabulated  together  with  some 
other  results  in  Tables  3.1  and  3.2. 

In  the  numerical  computation  the  slope,  da/de,  of  the 
true  stress-strain  curve  will  be  needed.   Therefore  the  value 
of  dS/de  will  be  converted  to  da/de  by  the  relation: 


da/de  =  (1  +  e )  "^  dS/de  +  (1  +  e)S 


3.1 
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3 . 2   Incremental  Wave  Test  Results-Oscilloscope  Photographs 


The  incremental  wave  test  for  each  specimen  contains 
about  20  runs  at  increasing  prestrain  levels  during  the 
quasistatic  loading  in  plastic  deformation.   All  the  velocity 
profile  traces  on  the  oscilloscope  screen  were  recorded  by 
Polaroid  pictures.   In  this  section,  some  typical  examples 
of  oscilloscope  trace  photographs  will  be  shown.   Since  the 
amount  of  data  for  7  aluminum  and  4  copper  specimens  is 
around  170  pictures,  data  reduction  was  done  by  using  a  com- 
puter program,  which  will  be  described  in  the  next  section. 
The  dependence  of  the  incremental  wave  behavior  on  both  pre- 
strain and  prestrain  rate  will  be  discussed  in  Section  3.4. 

Figures  3.3  and  3.4  show  examples  of  oscilloscope 
traces  of  incremental  particle  velocity  for  an  aluminum  speC' 
imen.  It  can  be  seen  in  the  case  of  small  prestrain  values, 
like  those  in  Figure  3.3(a)  and  (b),  that  the  velocity  pro- 
file shows  an  initial  step,  which  propagates  with  very    lit- 
tle attenuation  and  is  followed  by  a  more  slowly  rising  part 
that  attenuates  rapidly  along  the  gage  stations.   For  the 
cases  with  prestrain  larger  than  0.006,  the  rising  part  was 
flattened  out  even  in  the  trace  for  the  first  gage  station. 

Figures  3.5  and  3.6  show  examples  of  oscilloscope 
traces  for  incremental  particle  velocity  in  copper  specimens 
In  these  pictures,  it  can  be  seen  that  the  shape  of  the  ve- 
locity profile  varies   gradually  from  concave  upward  and 
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Pres trai  n 


e   =  0.0015 

0 


e   =  0.0060 

0 


e.  =  0 .0100 


Vertical  scale:  0.5  mv/cm,   Time  scale  (Nomi  nal  ly )  :  50  ;:sec/cm 
Figure  3.3   Oscilloscope  Trace  Records  (Aluminum  Specimen  No. 4) 
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Prestrain 
e    =  0.0200 


e     =0.0342 


e     =   0.0460 

0 


Vertical   Scale:   0.5  mv/cm.     Time  Scale   (Nominany):   50  ysec/cm 
Figure    3.4      Oscilloscope   Trace    Records  (Aluminum   Specimen    No. 4) 
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Prestrain      (aj    e      =    0.0009  (b)      e      =   0.0030 


(c)   e^  =  0.0070 


(d;   e^  =  0.0180 


Vertical  Scale  :  Station  No.  1  is  0.5  mv/cm 

Others  are  0.2  mv/cm 
Time  Scale  :  (Nominally)   100  vsec/cm 


Figure  3.5   Oscilloscope  Trace  Records  (Copper  Specimen  No.  4) 
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Prestrain   (a)  e   =  0.0220 

0 


(b)   e^  =  0.0330 
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(c)   e^  =  0.0400 


(d)   e„  =  0.04! 


Vertical  Scale  :  Station  No.  1  is  0.5  mv/cm 
Others  are  0.2  mv/cm 
Time  Scale  :  (Nominally)   100  psec/cm 


Figure  3.6   Oscilloscope  Trace  Records  (Copper  Specimen  No. 4] 
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attenuating  fast  to  concave  downward  and  attenuating  slower 
as  the  value  of  prestrain  increases. 

From  a  rough  examination  it  was  found  that  this  kind 
of  shape  changing  in  the  velocity  profile  will  give  a  family 
of  wave-propagation  speed  functions.   In  order  to  get  a 
quantitative  result,  the  wave  speed  functions  for  all  the 
200  incremental  tests  were  obtained  numerically  by  the  com- 
puter program. 

3 . 3   Velocity  Data  Reduction  and  Have  Speed  Function 


A  Fortran  program  was  written  to  reduce  the  large 
amount  of  incremental  wave  test  results  and  obtain  wave  speed 
functions  by  numerical  computations.   The  program  is  given  in 
Appendix  A-  1  .     Only  some  basic  procedures  are  explained 
briefly  in  the  following. 

(a)   Data  reading  and  conversion  into  velocity  units. 
Coordinates  of  about  15  points  of  each  velocity 
trace  on  the  oscilloscope  photograph  were  read 
and  converted  into  values  in  velocity  units  by 
using  the  velocity  transducer  calibration  factors 
obtained  from  a  subroutine  [see  Section  2.4]. 
( b  )   I nterpol ati  on . 

By  applying  an  i nterval -by- i nterval  cubic  Lagrange 
interpolation  subroutine,  the  values  of  velocity 
at  much  smaller  time  intervals,  At,  were  obtained. 
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At  was  chosen  to  be  4  usee  for  aluminum  results 

and  6  usee  for  copper  results. 

(c)   Computing  wave-propagation  speed  function. 

From  the  result  of  (b),  each  velocity  profile  is 

actually  composed  of  discrete  points  each  with  At 

time  interval.   These  data  points  are  expressed  by 

two  corresponding  arrays,  say  V(I)  and  T(I).   To 

obtain  wave-propagation  speed  of  velocity  level 

v,  =  k  Av,  a  routine  will  locate  the  interval  in 
k 

which  V,  lies,  and  then  use  linear  interpolation 
between  the  two  adjacent  points  to  obtain  the  cor- 
responding time  t.  and  store  it  in  another  array 
Q(k).   Wave-propagation  speed  computed  between 
gage  1  and  2  can  then  be  obtained  from 


c(v^)  =  distance/[Q2(k)  -  Q^  (k)] 


(3.2) 


In  this  manner,  wave-propagation  speeds  were  cal- 
culated between  velocity  stations  1  and  2,  2  and 
3,  3  and  4,  and  1  and  4.   Av  used  in  the  computa- 
tion is  equal  to  0.4  cm /sec. 
(d)   The  results  of  particle  velocity  and  wave- 
propagation  speed  are  printed.   They  are  plotted 
to  give  a  better  impression. 
Some  typical  examples  of  resulting  vel oci ty- time  pro- 
files of  an  incremental  wave  plotted  by  computer  are  given 
in  Figures  3.7  to  3.21.   In  these  figures  are  shown  also 
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results  from  the  theoretical  computation  of  Chapter  4. 

Figures  3.7  to  3.13  are  selected  from  specimen  No.  4 
of  copper.   Figures  3.14  to  3.21  are  selected  from  the  re- 
sults of  aluminum  specimen  No.  4.   Most  of  the  original  os- 
cilloscope traces  for  these  results  have  already  been  shown 
in  Figures  3.3  to  3.6. 

It  can  be  seen  from  these  figures  that  most  of  the 
results  were  good  except  that  the  gage  station  No.  2  of  the 
aluminum  case  shows  a  discrepancy  consistently.   In  general 
the  discrepancy  may  be  caused  by  the  following  possible 
sources  of  error:   (1)  error  in  reading  traces  on  the  photo- 
graph, (2)  error  in  magnetic  field  mapping,  (3)  error  in  lo- 
cating the  transducer's  position  in  the  magnetic  field  dur- 
ing the  wave  test,  (4)  error  in  the  calibration  of  the  os- 
cilloscope vertical  scale,  (5)  noise  signal  caused  by  the 
random  movement  of  wire  in  the  velocity  transducer's  hole, 
which  was  drilled  a  little  too  large.   For  the  case  just 
mentioned  it  was  No.  (5)  that  caused  the  error. 

The  elastic  bar-wave  propagation  speed  c   was  ob- 
tained by  measuring  the  time  interval  between  leading  edges 
of  incremental  velocity  profiles.   The  resolution  in  deter- 
mining the  position  of  the  leading  edge  is  very    poor  in  gen- 
eral.  It  can  be  seen  from  Figures  3.3  to  3.6  that  the  dis- 
tance between  two  adjacent  leading  edges  on  the  picture  is 
around  5  mm.   A  0.5-mm  error  will  produce  about  10%  error 
in  wave-propagation  speed.   No  systematic  variation  of  the 
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value  of  c   among  different  prestrain  levels  was  found.   The 

0 

result  of  c   is  the  averaged  values  of  76  cases  for  copper 
and  100  cases  for  aluminum.   Young's  modulus  was  calculated 
from  E  =  PC. 

5  12        2 

Aluminum:   c   =  5.13x10   cm/sec,   E  =  0.71x10    dyne/cm 

5  12        2 

Copper:   c   =  3.67x10   cm/sec,   E  =  1.20x10    dyne/cm 


By  checking  all  the  computer  output  of  wave-propagation 
speed  functions,  it  was  noted  that  the  variations  among  the 
curves  computed  from  gage  stations  1-2,  2-3,  3-4  and  1-4 
were  insignificant  and  random.   Therefore  it  was  concluded 
that  the  propagation  speed  of  a  given  level  of  incremental 
velocity  was  independent  of  propagation  distance.   The  wave 
speed  function  computed  from  gage  stations  1-2  hence  was 
used  to  express  the  wave  speed  function. 

From  the  computer  outputs  it  was  also  noticed  that 
wave-propagation  speed  functions  were  different  for  incre- 
mental waves  generated  at  different  prestrain  levels.   Fig- 
ures 3.22  and  3.23  show  the  typical  examples  of  these  re- 
sults for  copper  and  aluminum,  respectively.   It  was  appar- 
ent from  these  figures  that  the  value  of  wave  speed  at  a 
certain  velocity  level  is  usually  larger  for  larger  pre- 
strain values.   But  there  are  some  exceptions  like  curve 
No.  5  in  Figure  3.22  and  No.  7  in  Figure  3.23.   These  are 
not  the  only  exceptions,  since  many  curves  have  not  been 
plotted  in  the  figure  to  avoid  curves  mingling  together. 
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These  phenomena  were  examined  more  fully  by  correlating  para 
meters  obtained  from  wave  speed  curve  fitting  with  both  pre- 
strain  and  prestrain  rate  values. 


3 . 4   Wave-Propagation  Speed  Curve  Fitting  and  the  Discovery 
of  Prestrain  and  Prestrain  Rate  Dependence  of  the 
Incremental  Waves 


In  order  to  use  some  simple  parameters  to  character- 
ize the  family  of  wave-propagation  speed  functions,  various 
kinds  of  curve  fitting  techniques  were  tried,  among  which 
the  exponential  function  of  equation  (3.3)  below  was  found 
to  be  the  most  adequate  one.   Although  the  two  families  of 
c(v)  curves  for  copper  and  aluminum  look  so  different  in 
Figures  3.22  and  3.23,  they  could  both  be  fitted  well  by 
this  same  exponential  function  ,  with  different  parameter 
values.   It  was  found  from  the  computer  outputs  that  the 
wave-propagation  speed  lies  between  c  =  c   and  c  =  c  ,  where 
c   =  /T/ p  a  n  d  c   =  /Td  o / d  e ) / p  .   (In  finite  wave-propagation 


c   is  defined  as  c   =  /( dS/de )/ p .   For  the  incremental  wave 
the  initial  dimensions  corresponding  to  true  stress  and 
strain  values  are  the  initial  conditions.)   The  value  of  c 
decreases  as  the  prestrain  value  is  increased.   The  expres- 


sion c  =  (c  -  c  )/(c 
^      p     0 


c  )  was  used  as  the  nondimensional 
P 


wave-propagation  speed.   c  varies  from  1  to  0  as  c  varies 

from  c   to  c  .   The  incremental  velocity  was  also  made  non- 
0     p  -^ 

dimensional  by  dividing  it  by  a  velocity  parameter  v  .   The 
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form  of  the  curve-fitting  function  was 


c  =  exp[-(v/v^)^'] 


(3.3) 


in  which  v   and  n  were  determined  from  experimental  data. 
This  was  done  by  taking  the  natural  logarithm  twice  on  both 
sides  of  equation  (3.3),  which  was  then  reduced  to  a  linear 
equation  in  variables  ln[ln(l/c^)]  and  ln(v).   The  coeffi- 
cients of  this  linear  equation  were  then  determined  by  lin- 
ear regression.   A  simple  computer  program  was  written  to 
do  the  computations  for  determining  the  values  of  v   and 
n  for  all  the  170  wave-propagation  speed  curves. 

Some  typical  examples  of  the  results  of  the  curve 
fitting  are  shown  in  Figures  3.24  and  3.25.   It  can  be  seen 
that  the  fitted  curves  match  very    well  with  the  experimental 
data.   This  kind  of  good  quality  of  curve  fitting  was  also 
true  for  all  the  other  copper  data  and  65%  of  the  other  alu- 
minum data.   Tables  3.1  and  3.2  give  the  resulting  values 
of  V   and  n  together  with  corresponding  prestrain  rate  e 
and  plastic  wave  speed  for  copper  specimen  No.  4  and  alumi- 
num specimen  No.  4,  respectively. 

The  parameter  v   expresses  the  degree  of  outwardness 
of  the  c  curve  from  the  origin.   The  larger  the  value  of  v 
is,  the  farther  out  the  c  curve  lies.   Parameter  n  expresses 
the  curvature  of  the  c  curve.   Since  the  value  of  n  is  ^^ery 
sensitive,  the  resulting  value  of  n  from  curve  fitting  is 
not  as  reliable  as  is  the  value  of  v.   This  can  be  shown 
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by  taking  the  partial  derivative  of  equation  (3.3)  with  re^ 
spect  to  n  .   We  obtai  n 


Ac/c  -  -[ln(v/v^)](v/v^)"  An 


(3.4) 


For  example,  if  we  choose  n  =  1.0  for  copper  and  v/v   =  0.8, 

a  10%  variance  in  n,  i.e.,  An  =  0.1,  will  give  Ac/c^  =  1.8%. 

For  aluminum  the  value  of  n  is  around  5.   If  v/v   =  0.8 

again  and  An  =  0.5,  then  Ac/c  =  3.6%. 

The  values  of  v^  and  n  versus  the  prestrain  value  e 
c  0 

for  copper  specimen  No.  4  are  plotted  in  Figures  3.26  and 
3.27,  respectively.   At  first  glance  the  zigzag  relation  be- 
tween v^  and  prestrain  e   seems  to  mean  that  the  results  of 
wave-propagation  speed  contain  some  errors.   By  examining 
the  relation  carefully  and  correlating  it  with  valuesof  pre- 
strain rate,  it  was  found  that  whenever  the  value  of  e   in- 

0 

creased  the  value  of  v   dropped.   This  phenomenon  was  checked 
on  the  results  of  three  other  copper  specimens  and  five 
other  aluminum  specimens,  which  also  showed  the  same  char- 
acter.  Thus  it  was  concluded  that  the  character  of  small 
incremental  waves  propagated  on  the  quas i s tati cal ly  deformed 
rod  in  the  plastic  region  is  influenced  by  both  prestrain 
and  prestrain  rate . 

At  this  stage,  a  quantitative  relation  among  the  para- 
meter v^  and  prestrain  e   and  prestrain  rate  e   was  sought. 
It  was  found  that  on  the  e   and  loq(e  )  plane,  shown  in  Fiq- 

0  0 

ure  3.28,  the  contour  lines  of  constant  value  of  v   are 

c 
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close  to  a  family  of  parallel  straight  linrjs.   The  follow- 
ing linear  relation  was  assumed 


"^^  -  ^°9l0(%/^S)  =  h  ^  +  ^2 


(3.5 


where  m  is  the  slope  of  the  contour  lines  in  Figure  3.28, 

-  5    - 1 
e*  =  10    sec    is  an  arbi  trari  ly  chosen  reference  value  of 

the  prestrain  rate,  and  k-,  and  k^  are  constants.   Values  of 

m,  k-,  and  k^  were  determined  by  the  least  square  method. 

The  left-hand  side  of  equation  (3.5)  was  defined  as  a 

precondition  parameter  E,    =    me      -    loq(e  /e*)  in  an  attempt 

0  0   0 

to  express  the  dependence  on  e   and  e   by  a  single  parameter 
The  results  for  five  aluminum  and  four  copper  specimens  are 
listed  in  Table  3.3.   Results  of  copper  specimen  No.  4  are 
also  plotted  in  Figure  3.29. 

It  was  found  that  the  single  parameter  correlates  the 
results  for  any  one  specimen  fairly  well,  as  in  Figure  3.29. 
But,  as  can  be  seen  from  Table  3.3  the  value  of  m  in  the 
precondition  parameter  definition  was  not  constant  for  all 
specimens  of  one  material. 

It  was  not  found  possible  to  determine  any  correla- 
tion of  the  parameter  n  with  the  corresponding  prestrain  and 
prestrain  rate  values. 
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The  decimal  point  of  each  value 
of  V  (cm /sec)  indicating  the 
location  on  the  e^  and  logle^) 
p  1  ane . 


>  I  I  I...   I,  .1 I — I — I — I — 

0.00      0.01       0.02 


■  I  I 1 I I — I — 1 — I — I — I — 1- 


0.03 
Pres  tra i  n  e 


0.04 


Figure  3.28   Values  of  v^  on  the  e^  and  logle^)  PI 
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The  numbers  besides  each  point 
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Figure  3.29   Relation  Between  v   and  the  Precondition 
Parameter  E, . 
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Table    3.1       Values    of   Prestrain    Rate    and   Wave    Speed 

Function    Parameters    (Copper    Specimen    Ho. 4) 


^0 

eo' 
(D* 

4.02 

10"^/ 
(2)t 

3.51 

sec 

Ave. 

3.8 

cm/sec 
1  .53 

n 

5  'P 
10   cm/sec 

0.0009 

0. 

707 

0.7887 

0.0020 

2.38 

2.12 

2.3 

2.04 

0. 

.663 

0.6081 

0.0030 

1.12 

1  .07 

1  .1 

2.58 

0. 

.688 

0.5803 

0.0040 

3.93 

3.89 

3.9 

2.07 

0. 

.597 

0.5553 

0.0050 

2.55 

2.50 

2.5 

2.73 

0. 

.606 

0.5374 

0.0060 

1.37 

1.48 

1  .4 

3.75 

0. 

.  709 

0.5257 

0.0070 

-- 

1  .82 

1  .8 

3.85 

0. 

,624 

0.5164 

0.0080 

2.59 

3.09 

2.9 

3.71 

0. 

.600 

0.5096 

0.0090 

1  .00 

1  .35 

1  .2 

4.98 

0. 

.686 

0.5031 

0.0100 

3.67 

5.  15 

4.4 

4.88 

0, 

.619 

0.4979 

0.0121 

7.90 

8.70 

8.3 

3.18 

0, 

,683 

0.4896 

0.0140 

1.17 

2.03 

1  .6 

5.03 

0, 

.847 

0.4824 

0.0160 

4.85 

4.57 

4.7 

4.30 

0, 

.794 

0.4741 

0.0180 

2.42 

3.  12 

2.8 

5.47 

0. 

,914 

0.4665 

0.0200 

1  .50 

2.04 

2.0 

5.90 

0. 

,923 

0.4593 

0.0220 

3.00 

2.24 

2.5 

6.03 

0. 

,923 

0.4530 

0.0240 

2.92 

2.31 

2.6 

6.69 

0, 

,971 

0.4473 

0.0270 

5.  53 

4.87 

5.2 

6.44 

0, 

,925 

0.4368 

0.0290 

3.42 

4.  18 

3.8 

7.01 

0. 

.946 

0.4338 

0.0310 

3.42 

3.74 

3.6 

7.83 

1  , 

.007 

0.4283 

0.0330 

2.  25 

2.97 

2.6 

8.78 

1  , 

.041 

0.4228 

0.0350 

3.93 

3.98 

4.0 

8.35 

0. 

.988 

0.4182 

0.0376 

5.27 

5.53 

5.4 

8.16 

0, 

.966 

0.4110 

0.0400 

3.09 

4.41 

3.8 

9.43 

1  . 

,041 

0.4064 

0.0420 

0.42 

1  .02 

0.7 

10.59 

1  . 

,060 

0.4033 

0.0434 

0.42 

1  .16 

0.8 

14.27 

1  . 

,  358 

0.3990 

0.0488 

0.08 

0.03 

0.05 

18.71 

1  , 

,  157 

0.3904 

From    strain-time    records. 
From    load-time    records. 
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Table    3.2      Values    of    Prestrain    Rate    and    Wave    Speed 
Function    Parameters    (Aluminum  Specimen  No. 4) 


e 

0 

(D* 

.  10-^ 

(2)'^ 
2.6 

/sec 
Ave. 

2.6 

Vc 
cm/sec 

10.1 

n 
6.2 

10  ^  cm/sec 

0.0015 

1  .0650 

0.0030 

-- 

2.4 

2.4 

9.7 

7.8 

0.9345 

0.0060 

17.2 

16.0 

16.6 

11.4 

5.1 

0.8087 

0.0080 

5.2 

6.2 

5.7 

11.4 

5.8 

0.7522 

0.0100 

4.0 

5.2 

4.6 

12.6 

6.4 

0.7028 

0.0120 

3.6 

4.4 

4.0 

14.7 

5.5 

0.6606 

0.0140 

2.4 

4.1 

3.3 

14.3 

5.3 

0.6302 

0.0160 

9.8 

5.1 

7.5 

14.6 

5.6 

0.6013 

0.0200 

-- 

0.2 

0.2 

17.5 

4.6 

0.5465 

0.0250 

4.8 

5.7 

5.3 

16.9 

7.4 

0.5023 

0.0280 

3.7 

4.1 

3.  9 

21  .7 

4.8 

0.4711 

0.0300 

1  .6 

1  .8 

1  .7 

22.8 

4.0 

0.4563 

0.0342 

-- 

14.2 

14.2 

16.4 

5.2 

0.4243 

0.0400 

-- 

5.9 

5.9 

23.8 

3.6 

0.3835 

0.0460 

3.5 

3.5 

24.  1 

6.8 

0.3163 

From    Strain-time    records 
'''From    load-time    records. 


Table  3.3   Parameters  for  Functional  Representation 

Giving  v   in  Terms  of  e„  and  e  . 
^   C  0        0 


Specimen 

No. 

m 

K-,  (sec/cm) 
0.11 

-2.63 

Data  Numbers 

Aluminum  - 

-  1 

20.8 

6 

3 

11.4 

0.34 

-0.91 

14 

4 

48.4 

0.16 

-2.15 

15 

6 

85.8 

0.25 

-3.39 

13 

7 

34.2 

0.08 

-1  .42 

12 

Copper 

-  1 

96.3 

0.44 

-1  .58 

6 

2 

43.2 

0.21 

-0.72 

22 

3 

114.8 

0.61 

-1  .34 

21 

4 

46.4 

0.23 

-0.85 

27 

CHAPTER   4 
INTERPRETATION  OF  EXPERIMENTAL  RESULTS  BY  VARIOUS  THEORIES 


4 . 1   Equation  of  Motion  and  Compatibility  Condition 
for  Incremental  VJaves  of  Small  Amplitude 


The  governing  equations  for  one-dimensional  longitu- 
dinal wave  propagation  ina  bar  are  derived  under  two  general 
assumptions:   (1)  plane  cross  sections  remain  plane  and  (2) 
lateral  inertia  effects  are  negligible.   It  appears  that 
both  assumptions  are  fulfilled  in  the  present  experimental 
arrangement  for  the  incremental  wave  considered.   That  the 
plane  cross  sections  remain  plane  can  be  verified  directly 
from  the  experimental  results  for  the  elastic  wave  shown  in 
Figure  2.5  of  Section  2.3.   The  strain-gage  output  indicates 
response  from  the  surface  while  the  velocity  transducer  in- 
dicates  average  response  across  the  diameter  of  the  cross 
section.   The  identity  in  these  two  signals  verifies  that 
the  longitudinal  particle  velocity  is  uniform  over  a  given 
transverse  plane.   Lateral  inertia  effects  have  been  shown 
to  be  negligible  after  a  wave  has  propagated  a  distance  of 
one  diameter  from  the  impact  location  [see  Bell  (1966),  De 
Vault  (1965)].   In  the  present  case,  the  first  gage  station 
was  located  15  cm  or  15  diameters  from  the  impact  collar  at- 
tached on  the  rod.   This  first  gage's  output  was  used  as  the 
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boundary  condition.   Therefore  the  effects  of  lateral  iner- 
tia are  negligible  in  the  region  following  the  first  gage 
station. 

Lagrangian  coordinates  based  on  the  initial  dimen- 
sions of  each  prestrained  configuration  are  used  for  de- 
scribing the  small  incremental  waves  generated  in  that  pre- 
strained configuration.   The  coordinate  x  =  0  is  chosen  as 
the  position  of  the  first  gage  station  at  the  instant  when 
the  leading  edge  of  the  incremental  wave  reaches  the  first 
gage  station.   The  distance  between  adjacent  gages  at  this 
time  (t  =  0)  is  equal  to  12.7(1  +  e  )  cm,  where  e   is  the 
value  of  nominal  prestrain.   Let  u(x,t)  be  the  displacement 
at  time  t  of  the  cross  section  that  was  at  position  x  at 
t  =  0.   Then 


ju 
3x 


(4.1) 


9t 


(4.2) 


where  e  is  the  incremental  strain  and  v  is  the  incremental 
particle  velocity  at  the  section  under  consideration.   In 
this  chapter  the  symbols  a,  e  and  v  will  be  used  to  repre- 
sent the  incremental  quantities  of  stress,  strain  and  ve- 
locity.  Tensile  incremental  stress  and  strain  are  reckoned 
positive.   Displacement  and  particle  velocity  in  the  posi- 
tive x-direction  are  reckoned  positive.   Therefore  in  the 
present  case,  a  tensile  incremental  wave  is  accompanied  by 
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negative  particle  velocity.   From  equations  (4.1),  (4.2) 
the  compatibility  condition  can  be  obtained  as 


de 
dt 


(4.3) 


The  equation  of  motion  for  an  element  of  the  bar  can  be  de- 
rived from  the  impulse-momentum  principle  as 


9o     8v 
3x    ^  dt 


(4.4 


The  material  density  p  is  assumed  to  be  constant  and  inde- 
pendent of  the  prestrain  value. 

4 . 2   Boundary  and  Initial  Conditions 


The  boundary  condition  for  the  theoretical  solutions 
considered  in  this  chapter  was  taken  from  the  experimental 
particle  velocity  data  of  the  first  gage  station.   Accord- 
ing to  the  sign  convention  mentioned  in  Section  4.1  the 
value  of  the  velocity  should  be  negative  in  the  theoretical 
computations. 

The  initial  conditions  for  the  incremental  stress, 
strain  and  velocity  are 


0 


for   t  -  0,  X  >  0     (4.5) 


These  initial  conditions  are  appropriate  since  in  the 

-4 
duration  of  10    sec  (more  than  the  time  needed  for  the 
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incremental  wave's  leading  edge  to  travel  through  the  four 

gage  stations)  the  prestress  and  prestrain  of  the  quasistatic 

-5    - 1 
loading  at  a  strain  rate  of  the  order  of  10    sec    can  be 

considered  constant  and  the  initial  particle  velocity  can 

be  consi  dered  zero . 


4 . 3   Stra i n-Rate- I ndependent  Theory 

From  the  experimental  results  of  Section  3.3,  the 
wave-propagation  speed  for  a  given  level  of  particle  velo- 
city is  a  constant  value  along  the  propagation  distance 
through  the  four  gage  stations,  as  would  be  predicted  by 
the  simple  wave  solution  of  the  one-dimensional  strain-rate- 
independent  theory.   An  exception  occurs  in  the  vicinity  of 
the  plateau  in  the  particle  velocity  records  for  aluminum. 
Since  the  curves  do  not  approach  the  same  particle  velocity 
plateau,  it  is  apparent  that  no  strain  plateau  is  approached 
in  the  time  of  the  experimental  record.   A  very    complete  ac- 
count of  the  rate- independent  theory  for  finite   amplitude 
plastic  wave  propagation  is  included  in  Cristescu  (1967). 
In  this  section,  the  basic  analytical  method  of  the  strain- 
rate-independent  theory  will  be  applied  similarly  for  incre- 
mental waves  by  assuming  that  there  exists  a  unique  dynamic 
incremental  s tress -strai n  relation 


o  -  h{e) 


(4.6) 
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for  a  given  incremental  wave  at  a  certain  prestrain  value. 
The  dynamic  response  function  is  further  assumed  to  be 

single-valued  and  to  satisfy  the  conditions  dh/de  >  0  and 

2     2 
d  h/de  <_   0.   The  present  problem  is  of  an  inverse  type, 

i.e.,  to  find  the  function  h{c)  from  the  experimental  data 

of  the  particle  velocity. 

Introducing  equation  (4.6)  into  (4.4)  transforms  the 

equation  of  motion  to 


1  dji  iG 
p  de  9x 


(4.7) 


Equation  (4.7)  and  the  compatibility  equation  (4.3)  form  a 
system  of  first  order  quasilinear  partial  differential  equa 
tions.   The  characteristics  are  given  by 


d_x 

dt 


c(e) 


(4.8) 


where 


p  de-" 


4.9 


The  differential  relations  satisfied  along  the  characteris 
tics  are 


dv  =  ±  c ( e)  de 


(4.10) 


The  upper  and  lower  signs  in  equations  (4.8)  and  (4.10)  cor- 
respond to  each  other.   It  can  be  shown  that  the  initial 
constant  state  of  initial  condition  (4.5)  gives  a  simple 
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wave  solution.   The  characteristics  of  positive  slope  are  a 
family  of  straight  lines,  and  a,  e,  and  v  are  constant  val- 
ues along  each  positive  characteristic.   The  differential 
relation  (4.10)  becomes 


dv 


-  c ( e)  dc 


(4.11) 


where  the  differentiation  can  be  taken  along  any  direction 
in  the  simple  wave  domain. 

The  incremental  wave-propagation  speed  as  a  function 
of  particle  velocity  was  obtained  from  experimental  data. 
The  incremental  response  function  h(e)  can  then  be  obtained 
from  it,  as  follows.   Combining  (4.9)  and  (4.11)  yields 


da 


PC  dv 


4.12 


Integrating  (4.11)  and  (4.12),  we  obtain 


e(  V  . 


/   dv/c ( V 

0 


(4.13 


and 


a(v.) 


-  /   PC ( V  )  dv 

0 


4.  14 


If  the  integration  in  both  equations  (4.13)  and  (4.14)  can 
be  carried  out  explicity,  then  a  closed  form  relation  be- 
tween a  and  e  is  obtained.   Otherwise,  the  one-to-one  cor- 
respondence between  a  and  e  can  be  obtained  by  integrating 


(4.13)  and  (4.14)  numerically  with  v.  increased  by  incre- 
ment Av  at  each  step. 

The  fitted  equation  of  wave-propagation  speed,  equa- 
tion (3.3),  can  be  rewritten  as 


c(v)  "  c   +  (Cq  -  Cp)  exp[-(v/v^)"] 


(4.15) 


where  values  of  c,  v  and  n,  constants  at  each  prestrain, 
were  given  in  Tables  3.1  and  3.2  for  copper  specimen  No.  4 
and  aluminum  specimen  No.  4,  respectively. 

The  integration  was  carried  out  by  using  a  recursion 
formula.   From  equations  (4.13)  and  (4.14),  the  numerical 
schemes  are  simply 


S-  +  1  =  "(%-  +  l)  -     S-  +  Av/c(v".^^) 

i  =  1  ,  2 ,  ... 


(4.16) 
(4.17) 


where  Av  = 


i  +  1 


i  +  1 


2^%-  +  %-  +  l)'  ^1 


0, 


v.  =  (i-1)  Av.  (flote  that  the  value  of  v  is  negative, 


de  =  -  dv/c  =  d|v|/c,  etc.). 

In  practice  Av  =  0.1  cm/sec  was  used  and  the  compu- 
tation was  carried  out  to  v  =  40  cm /sec.   Convergence  was 
checked  by  comparing  with  results  obtained  from  Av  =  0.2  cm/sec: 
differences  in  the  computed  values  of  a  and  e  at  v  =  40  cm/sec 
were  less  than  0.001%.   Figures  4.1  and  4.2  show  some  exam- 
ples of  the  incremental  stress-strain  curves  at  various 
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prestrain  levels  for  copper  specimen  No.  4  and  aluminum 
specimen  No.  4,  respectively. 

4 . 4   Strain-Rate-Dependent  Theory 

Since  the  first  use  of  a  semilinear  rate-type  consti- 
tutive equation  by  Sokolovsky  (1948)  and  Malvern  (1949)  and 
the  introduction  of  a  quasi  linear  equation  by  Malvern  (1949), 
several  other  kinds  of  rate- type  constitutive  equations  have 
been  investigated  in  many  published  papers  [see  Cristescu 
(1963),  (1972),  Lubliner  (1964),  Lubliner  and  Valathur  (1969)], 

In  this  section,  two  rate-type  constitutive  equations 
in  terms  of  incremental  stress  and  strain  were  considered 
and  numerical  solutions  were  obtained  for  incremental  waves 
at  several  prestrain  levels  of  one  aluminum  case  and  one 
copper  case.   The  results  show  that  the  simple  wave  solution 
can  also  be  approximated  in  certain  types  of  constitutive 
laws  of  the  strain-rate-dependent  theory. 

4.4.1   Constitutive  Laws  Considered 

The  general  form  of  the  quasilinear  constitutive 
equation  as  suggested  by  Cristescu  (1963)  is 


e    -    ^{o  ,z)a    +  ^\){o  ,z) 


(4.18) 


where  the  dot  represents  the  time  derivative  of  the  variable 


92 


at  a  material  point  x.   Two  special  forms  of  equation  (4.18) 
are  considered,  the  first  is  the  semilinear  model  of  Malvern^ 
type 


e  =  Fa  +  p-[a-Ae] 


(4.19) 


The  second  kind  is  a  quasilinear  model  proposed  by  the  pre- 
sent author  as 


ll   +  pe]i  +  |[a  -  Ae] 


4.20 


In  equations  (4.19)  and  (4.20),  the  parameters  p  and  k  are 
constants  to  be  determined,  E  the  Young's  modulus  and  A  the 
slope  of  the  quasi  static  true  prestress-prestrain  curve  at 
the  prestrain  level  where  the  incremental  wave  is  generated. 
The  term  Ae  represents  the  approximate  form  of  the  quasi- 
static  nominal  stress-strain  curve  based  on  the  initial  di- 
mensions of  the  prestrained  state  when  the  incremental  wave 
started.   (o  -  Ae)  is  the  overstress.   A  restriction  imposed 
on  the  instantaneous  term  ^{a,e)    was  suggested  by  Cristescu 
(1972).   It  was  that  whenever  a  <  0  or  a  <  Ae  the  value  of 
ct)(o,c)  should  be  switched  to  1/E. 


4.4.2   Characteristics  and  Associated  Interior  Differential 
Equations 

After  introducing  equation  (4.18)  into  equation  (4.3) 
and  then  combining  with  equation  (4.4),  we  can  write  the  re- 
sulting system  of  the  equations  in  the  following  matrix  form 


93 


0 


4) 


1 


0 


0    ro 


X  ) 


(4.21) 


-ij; 


where  the  subscripts  t  and  x  denote  partial  differentiation 
The  theory  for  systems  of  first  order  quasilinear  partial 
differential  equations,  as  represented  by  equation  (4.21), 
is  given  in  Jeffrey  and  Taniuti  (1964).   The  characteristic 
wave  speeds  associated  with  such  a  system  are  the  roots  of 
the  determinantal  equation 


PC 


1 


4.22) 


With  c  =  1  /  /p(^  the  characteristic  equations  are 


dx  =  ±  c  dt 


4.23 


The  null  vector  I      =  [i^i.  ^o]  for  equation  (4.21),  which 
corresponds  to  the  characteristic  wave  speed  c,  is  given  by 


[^1^  ^2^ 


PC 


[0,  0] 


(4.24) 


Choosing  an  arbitrary  scale  factor  such  that  I.    -    1  in  equi 
tion  (4.24),  we  get 


[-C(l),  1] 


(4.25) 
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The  interior  differential  equation  corresponding  to  the 
characteristic  given  by  dx  =  c  dt  is  given  by 

r 


[-c<t>,    1] 


P   0 
0    (b 


[-C(|),  1]  -(    >  dt    (4.26) 


which  is  rewritten  as 


do-pcdv=-Ydt, 


dx  ==  +  c  dt 


For  the  negative  characteristic  given  by  dx 
replace  (+c)  by  (-c)  in  equation  (4.27). 


(4.27) 


c  dt,  we 


4.4.3   Semilinear  Model  —  Numerical  scheme 

The  governing  differential  equations  given  in  (4.19) 
and  (4.27)  are  written  in  the  finite  difference  form  with 
4)  =  1/E  and  ij.'  =  k(o  -  Ae)/E  in  the  present  case.   The  char- 
acteristic wave  speed  c  of  equation  (4.22)  becomes  a  con- 
stant, equal  to  the  elastic  bar  wave  speed  c   =  /E/p.   The 
system  of  finite  difference  equations  is 


Ao 


E  Ae  =  -  [a  -  Aejk  At   along  dx  =  0 


4.28 


Aa  ±  pc   Av 

0 


[a  -  Ae]k  At   along  dx  =  ±  c^  dt  (4.29a,b 


Figure  4.3(i)  shows  the  solution  domain  in  the  x-t 
plane.   Mesh  point  types  involved  are  sketched  in  the  same 
figure  with  labels  (ii),  (iii),  (iv)  and  (v).   This  is  the 


95 


At 


r 


(i) 


Figure  4.3   Schematics  of  Solution  Domain  and  Mesh  Point 
Types  for  Semilinear  Model. 
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same  algorithm  as  that  given  by  Bianchi  (1963)  and  Banerjee 
(1972). 

The  unknowns  at  a  typical  point  are  stress  and  strain 
for  the  boundary  and  stress,  strain  and  velocity  for  a  gen- 
eral interior  point.   These  are  determined  by  solving  appro- 
priate selected  set  of  equations  from  equations  (4.28)  and 
(4.29a,b)  as  follows: 
General  boundary  point  [Ref.  Figure  4.3(iii)] 

Solution  is  known  at  points  b  and  c.   Solution  at 
point  d  is  obtained  from 

o,-o    +  pc  (v  ,  -  v  )  =  -  -J-k  [a,  +  o  -  A(e,  +  c  )]  At     (4.30) 
d    c    od    c     2'-dc     d    c         ^ 

Mesh  point  type  4.3(ii),  representing  the  first  point  to  be 
computed,  is  governed  by  equations  identical  to  those  above, 
except  that  At  in  equation  (4.30)  is  replaced  by  0.5  At. 

General  interior  point  [Ref.  Figure  4.3(v)] 

Solution  at  points  a,  b,  c  is  known.   Solution  at 
point  d  is  obtained  from 

a.-o     -  PC  (v  ,  -  v  )  =  -  |-k  [a,  +  a  -  A(c.  ,  +  e  )]  At     (4.32) 
d    aod    a'     2'-d    a    ^d    a 


0,-0     +pc(v,  -v)  =  -4-l<[o.  +  0  -A(e,  +  £)]At     (4.33 
d    co^d    c     2^d    c    ^d    c'-^ 


^d  -  %  -  ^(^d  -  %^   -  -  -^^^'d  '  %  -   ^^^d  "  %^^   ^^     (^-3^^ 
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Mesh  point  behind  leading  wave  [Ref.  Figure  4.3(iv)] 

The  difference  equations  are  the  same  as  equation 
(4.32)  through  (4.34),  except  that  At  in  equation  (4.33) 
is  repl aced  by  0.5  At . 

A  Fortran  IV  computer  program  was  written  to  perform 
the  computation.  Because  of  the  linear  prestress-prestrain 
function  assumed  in  the  overstress  term  of  equation  (4.19), 
direct  solution  from  the  linear  simultaneous  difference 
equations  is  possible.  In  order  to  make  the  actual  gage 
stations  fall  on  mesh  points  exactly  in  each  case  for  vari- 
ous prestrain  levels,  the  mesh  size  is  controlled  by  an  in- 
teger parameter,  denoted  by  NX,  which  represents  the  number 
of  intervals  between  two  adjacent  gage  stations.  Mesh  size 
dimensions  Ax  and  At  are  obtained  from 


Ax 


(1  +  e  )  d 
I; 0 

NX 


(4.35 


At  = 


(4.36) 


where  d  is  the  distance  between  two  adjacent  gage  stations 

at  zero  prestrain  and  e  is  the  value  ofprestrain.   NX  =  16 
^  0  ^ 

was  chosen  in  the  computations,  and  thus  the  value  of  At 
varies  from  0.7736  psec  to  0.8123  psec  as  the  prestrain  in- 
creases from  zero  to  0.05  in  the  aluminum  case.  At  ~    1.1  usee 
in  the  copper  case.   Boundary  data  were  read  at  constant 
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time  intervals  from  the  experimental  results.   These  were 
then  interpolated  by  using  the  cubic  Lagrangian  step-by-step 
routine  to  get  values  at  the  mesh  points  with  intervals  of 
At  obtained  in  each  case.   Convergence  of  the  computation 
was  checked  by  seeing  that  halving  the  mesh  step  (doubling 
the  value  of  NX)  did  not  produce  any  observable  change  in 
the  solution. 

Results  of  the  numerical  computation  will  be  discus- 
sed in  detail  in  Section  4.4.5.   Some  remarks  are  given  here 
In  general  the  semi  linear  model  of  the  Malvern  type   does 
not  give  a  simple  wave  solution.   Some  authors  recommended 
the  use  of  a  nonlinear  function  of  overstress  to  get  a  bet- 
ter result.   However,  Banerjee  (1972)  has  obtained  solutions 
for  the  semilinear  model  with  both  linear  and  exponential 
overstress  and  for  the  quasilinear  model  with  linear  over- 
stress  in  a  torsional  incremental  wave  case.   By  checking 
the  results  carefully,  it  was  found  that  the  result  of  the 
exponential  overstress  model  made  no  substantial  improvement 
over  that  of  linear  overstress.   It  was  the  quasilinear 
model  which  gave  the  best  result  among  the  three  models 
considered.   Therefore  a  quasilinear  model  was  considered 
next . 
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4.4.4   Quasi  linear  Model  —  An  I ntegro-Di f f erenti a1  Approach 
and  Numerical  Scheme 


The  behavior  of  the  theoretical  solution  for  the 
quasi  linear  constitutive  model 


ll   +  pe]a  +  |[o  -  Ae] 


(4.20) 


in  the  incremental  wave  propagation  will  be  considered.   The 
parameters  p  and  k  are  two  constants  to  be  determined  by 
matching  the  theoretical  solutions  with  the  experimental  re- 
sults for  incremental  waves  at  various  prestrain  levels. 
( i )   Remarks  concerning  the  choice  of  the  function  ^{a  ,z) 
The  introduction  of  the  instantaneous  plastic  re- 
sponse term  was  first  made  by  Cristescu  (1963)  on  a  purely 
phenomenological  basis,  in  order  to  give  a  better  descrip- 
tion of  the  dynamic  plastic  response  of  a  material.   How- 
ever, Holt  et  al .  (1967)  performed  tests  with  a  series  of 

commercial  aluminum  alloys  and  found  that  s tress -stra i n 

-2     -1       3 
curves  at  the  higher  strain  rates  of  10    sec    to  10 

sec'   were  all  coincident.   This  means  that  for  representing 
the  behavior  of  that  particular  material  in  that  high  strain 
rate  range,  the  existence  of  a  single  dynamic  stress-strain 
curve  can  be  assumed.   Most  authors  have  determined  the  dy- 
namic response  function  cj)  by  assuming  a  simple  form  of  dy- 
namic stress-strain  curve, suchas  a  =  Be"  [see  Cristescu  (1972), 
Banerjee  (1972)].   This  always  resulted  in  a  complicated 
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form  of  the  function  ({>  itself.   The  solution  to  the  conse- 
quent system  of  differential  equations  then  required  a  com- 
puter time  of  the  order  of  minutes  in  a  third-generation  com- 
puter.  The  cost  problem  became  critical  in  the  present  work 
since  solutions  for  several  cases  at  various  pre  strain  lev- 
els of  a  series  of  incremental  wave  tests  were  to  be  inves- 
tigated.  For  this  reason,  the  simple  form  of  (}>  =  1/E  +  pe 
was  proposed  and  a  solution  procedure  was  found  such  that 
the  computing  time  needed  to  obtain  a  solution  was  reduced 
to  about  half  a  minute.   The  function  cf)  =  1/E  +  pe  gives  an 
instantaneous  response  curve  of  the  form 

1 


a  =  -  1 n(l  +  pEe) 

P  r    / 


4.37) 


This  curve  is  qualitatively  satisfactory,  since  it  is  a 

2     2 
s i ngl e- val ued  function  with  9o/9e  >  0,  9  o/9e   <  0. 

A  restriction  proposed  by  Cristescu  (1972),  on  the 
form  of  (})  in  equation  (4.18)  is  that  the  slope  of  the  instan- 
taneous response  curve  must  lie  between  the  slopes  of  the 
elastic  curve  and  of  the  relaxation  boundary  at  the  same 
strain.   For  the  present  case  this  requires 


A  < 


1  +  pEe 


<  E 


(4.38) 


which  was  always  satisfied 
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( 1 i )  Formulation  of  the  interior  differential  equations 
in  an  i ntegro-di f f erenti al  form 

It  will  be  shown  that  the  quasilinear  constitutive 
equation  (4.20)  may  be  solved  for  e(t)  in  terms  of  o(t)  and 
t.   Thus,  one  of  the  unknown  variables  is  el i mi nated, and  the 
resulting  system  of  differential  equations  can  be  solved  in 
a  more  efficient  way. 

Consider  equation  (4.20),  which  represents  the  con- 
stitutive law  assumed  for  a  material  point  x. 


e  ^  Lf  +  Pe]o  +  f[o  -  Ac] 


4.20) 


The  initial  condition  is 


0 


at 


t  =  0 


(4.39) 


(In  this  discussion  time  t  =  0  represents  tlie  moment  when 
the  leading  edge  of  the  incremental  wave  just  arrives  at  x. 
This  differs  from  the  initial  time  of  Section  4.2).   The 
solution  to  equation  (4.20)  without  the  second  terms  on  the 
right-hand  side  may  be  written  as 


£(t) 


^.weP^(^ 


(4.40) 


where  w  =  constant.   The  solution  to  the  complete  equation 
(4.20)  may  then  be  given  by  the  variation  of  parameter  tech 
nique  [Hildebrand  (1962),  Suliciu  et  al .  (1974)]  as 
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e{t)    =  -  ^+  w(t)  eP^(t) 


(4.41) 


where  now  w(t)  is  a  variable  to  be  determined.   Introducing 
the  initial  condition  (4.39)  into  equation  (4.41),  gives 


w(0) 


pE 


He  make  the  change  of  variable 


w(t)  =  ^  +  w^ (t)  , 


w^ (0)  =  0      (4.42) 


and  substitute  equation  (4.41)  into  equation  (4.20)  to  get 


•   ^  Ak 
^  ^  T  ^1 


|([a(t)  +  ^]  e 


p  o  ( t )    _^_ 
PE 


4.43) 


Solution  to  equation  (4.43)  may  be  obtained  by  apply- 
ing the  variation  of  parameter  method  once  again.   The  re- 
sult i  s 


w^(t 


--^t  t  -^t 

F  e  ^   /  q(s)  ds  -  -^(1  -  e  ^  ) 


(4.44) 


whe  re 


q(s)  =  [a(s)  +  ^]  exp[^s  -  pa(s)] 


(4.45) 


Substituting  equations  (4.44)  and  (4.42)  in  equation  (4.41) 
we  get  an  explicit  representation  for  e  in  terms  of  a  and  t 


£(t)  =  p^  [F{a(t),  t}  -  1] 


4.46 
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where    F{a(t),t}     is     a    functional    depending    on    the    (unknown) 
history    of    o(t),    given    by 


F{o(t),t}    =   pk  /  [a(s)   +  ^]   exp{p[a(t)   -    o(s)]   -  -^(t  -   s)}  ds  + 

0 


exp[pa(t)  -  -^t] 


(4.47 


In  order  to  incorporate  it  with  the  equations  of  mO' 
tion  and  compatibility,  equation  (4.46)  is  differentiated 
with  respect  to  t,  giving 


!_£  ^  F{a,t }  la    k 
It      E    at    E  ° 


F{o,t}] 


4.48 


where  for  brevity  F{a(t),t}  has  been  written  as  F{a,t}.   By 
comparing  equation  (4.48)  with  the  general  form  of  quasi  lin- 
ear constitutive  equation  (4.48),  the  values  of  the  functions 
(t>   and  ii    in  equation  (4.21)  at  any  time  can  be  expressed  in 
terms  of  the  value  of  F{a,t}  as 


=  F{o,t}/E 


4.49) 


|a  +  ^2  n  -  F{o,t}] 


(4.50 


By  applying  the  results  of  Section  4.4.2,  we  obtain 
the  characteristic  wave  speed  of  equation  (4.22)  in  terms 
of  F{a  ,  t)  as 
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c    =    l//pF{a,t}/E 


or 


c    =    Cq/ZfToTTT 


(4.51) 


where  c   is  the  elastic  bar  wave  speed.   The  interior  dif- 
ferential equations  associated  with  the  characteristic  equa- 
tions 


dx 
dt 


/F  { o ,  t  } 


(4.52) 


may  then  be  obtained  from  equation  (4.27)  as 


do  :r  PC  dv  =  k  [A  (1  -  p-fi^)  -  ^^^-i    dt       (4.53) 


Only  two  unknown  variables  a  and  v  are  involved  in 
the  integration  of  equation  (4.53).   The  strain  is  obtained 
from  equation  (4.46)  after  solution  for  a  has  been  obtained. 
Some  difficulties  are  observed  from  solving  equation  (4.53), 
since  the  values  of  F{a,t}  are  not  known  a  priori,  because 
a(t)  is  not  known  and  equation  (4.47)  for  F{a,t}  contains  an 
integration  from  s  =  0  to  s  =  t  with  the  integrand  a  function 
of  a(s).   A  recursion  formula  of  numerical  integration  was 
found  that  decouples  the  history  dependence  of  o(t)  in  F{a,t} 
as  f ol 1 ows 


F„  -    exp(pa   -  at  )  +  pk  exp[p(a   -  a   ,)  -  a  At]{F   ,  + 
n     ^'^n     n    ^    '^'-^^n    n-1  n-1 


At[A.  +  ho      +    o      J]  exp[-l~(a  At  -  p(a   -  o      J)]} 
1    2'n     n-l"^    KL2V        ^^n     n-1 


(4.54 
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where  a  =  Ak/E,  A-,  =  A/pE,  At  =  t   -  t  _  ■,  and  the  subscripts 

n  and  n-l  represent  the  function  values  at  t  =  t   and  t  - 

'^  n 

t   -,,  respectively.   All  the  (  n-l  ) -s  ubscri  pted  quantities 

are  known  from  the  previous  iteration.   Therefore,  equation 

(4.54)  makes  the  value  of  the  function  F  depend  explicitly 

on  o  ,  the  present  unknown  value  of  stress.   This  makes  the 
n      ^ 

iteration  procedure  of  solving  for  o  and  v  at  a  mesh  point 
become  possible.   The  numerical  scheme  is  described  in  the 
fol  lowing. 
( i  i  i  )   Numerical  scheme 

The  method  of  integration  consists  of  integrating 
equation  (4.53)  along  the  corresponding  characteristic 
curves,  defined  by  equation  (4.52).   These  two  families 
of  characteristic  curves  have  variable  slopes  in  the  load- 
ing domain.   However,  the  grid  used  in  integration  was  a  re- 
gular one  with  constant  slopes  ±  c.   Figure  4.4  shows  sche- 
matically the  solution  domain  and  types  of  grids.   The  stand' 
ard  Courant- I saacson-Rees  procedure  [see  Jeffrey  and  Taniuti 
(1954)]  was  applied  and'the  convergence  of  iterations  is  as- 
sured, si  nee  c  <  c  . 
—   0 

General  interior  point  [Ref.  Figure  4.4(iii)] 

The  solution  is  known  at  points  a,  b  and  c.   Points 
L  and  R  are  intersections  of  lines  ab  and  be  with  the  lines 
Ld  and  Rd  tangent  to  characteristic  curves  at  point  d.   Fi- 
nite difference  formulations  of  equation  (4.53)  along  the 
characteristics  Cd  and  Rd  may  be  written  as 
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<w—  X 


(i) 


Ax/ 2 


(ii) 


( i  i  i  ) 


Figure  4.4   Schematics  of  Solution  Domain  and  Mesh  Point 
Types  for  Quasilinear  Model. 


107 


0,-0 

d 


pc  -.  -, 


1  ^    1  "d  ^  "l 


d    2 


]  kz  At 


(4.55) 


/  r  ,  d     d    2 


(4.56) 


where  A-,  =  A/pE  and  F,  is  unknown  and  depends  on  o,  as  given 
by  equation  (4.54).   The  meaning  of  parameter  z  can  be 
perceived  from  Figure  4.4(iii),  and  it  can  be  shown  that 


/F. 


1  +  yF, 


(4.57) 


The  variables  a^,  a,  ,  v^  and  v,  in  equations  (4.55)  and  (4.56) 
can  be  replaced  by  the  following  expressions  from  linear 
interpolation, 


R 


a,  +  r  (a 
b      c 


°b> 


o,  +  r  (o 
b      a 


v,  +  r  (  V 
b      c 


\    ^    '^\ 


4.58) 


where    r  =  bR/bc  -  bL/ba  =  2/(1  VF.) 


(4.59) 


Substituting  equation  (4.58)  into  equation  (4.55)  and  (4.56), 

we  obtain  the  solutions  for  o,  and  v, 

d      d 
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o^    -    -^    [20^^(1  -  r)  +  r(a^  +  a^)]  +  s^  -h  rs3(v^  -  vj  (4.60) 


'd    2 


[2vj^(l  -  r)  +  r(Vg  +  v^)]  +  rs^io^    -    o^) 


(4.61) 


where 


2F  ,  -  kz  At 

Q 

2Fj  +  kz  At 


(4.62a) 


2A^{F^  -  1 )  kz  At 


2F.  +  kz  At 
d 


(4.62b) 


c  /FT 
0   d 


2F^  +  kz  At 


(4.62c) 


2F,  -  kz  At 
d 

4  c^/FT 
0   d 


(4.62d) 


The  iteration  proceeds  by  assuming  first  o,  -    o , ^    =  a,  . 

The  values  of  F,  and  also  z,  r,  s-,,  s^,  s^,  s,  are  evaluated 
next.   First  solutions  a,    ,  v,     can  then  be  obtained 

from  equations  (4.60)  and  (4.61).   The  process  is  repeated  by 


using  0 


Xo 
n-1) 


(1)  , 


1 


A) 


(0, 


or  in  general  a,  =  Ao 


(n) 


(1  -  \)    ctj''"  '  '  ,    until  the  prescribed  agreement  is  obtained 
between  the  successive  a  ,  values.  The  parameter  A  is  used 
here  to  improve  the  rate  of  convergence,  since  otherwise 
some  kind  of  slowly  convergent  oscillation  occurs.   In  the 
final  computation,  A  =  0.9  was  used,  and  it  resulted  in  an 
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accuracy  of  |  [a  , 


(n) 


(n-l)n/,  (n),  , 


-  a,^"  '''\/a.^"'\     <   0.00005  in  an  aver- 


age of  4  iterations.   If  A  =  1.0  is  used,  i.e.,  a,  =  a, 

d    d 


(n) 


the  number  of  iterations  required  is  around  25. 
General  boundary  points  [Ref.  Figure  4.4(ii)] 

The  only  unknown  is  the  value  of  stress  at  point  d. 
The  particle  velocity  values  are  given  on  the  boundary  points 
The  iteration  scheme  is  the  same  as  that  of  the  general  in- 
terior points,  except  that  equations  (4.60)  and  (4.61)  are 
replaced  by 

^d  "  ^^3^^\  '   ^d^  ^  ""^V  "  ^d^^  ^  ^l^^^^^  -  "^^  ^   '^'^c^  "^  ^2     (4.63) 

The  mesh-size  controlling  parameter  was  arranged  in 
a  way  similar  to  that  used  in  the  semilinear  model.   The 
values  of  At  used  were  around  2.0  ysec  for  both  copper  and 
aluminum  specimens.   Halving  of  this  step  size  did  not  pro- 
duce any  significant  change  in  results  but  increased  the 
computer  time.   The  total  computer  time  on  the  IBM  system 
360/70  was  around  30  seconds  for  a  loading  duration  of  350 
ysec.  Computer  program  is  given  in  Appendix  A-2. 

4.4.5   Numerical  Results 


( i )   Resul ts  for  copper 

Experimental  results  of  specimen  No.  4  were  used  for 
comparison  with  the  numerical  solutions  by  theoretical  mod- 
els.  It  was  found  that  some  substantial  discrepancies  exist 
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between  the  solution  by  the  semilinear  model  and  the  experi- 
mental data.   A  typical  example  is  given  in  Figure  4.5  for 
the  case  of  prestrain  e   =  0.0290.   The  discrepancies  are 
mainly:   (1)  the  v-t  curve  shapes  are  different  and  (2)  the 
propagation  speed  of  a  given  level  of  particle  velocity  is 
not  a  constant  according  to  the  numerical  solution  by  the 
semilinear  model  but  varies  by  as  much  as  about  30%  in  this 
case . 

Most  of  the  effort  was  therefore  spent  on  matching 
the  solutions  according  to  the  quasilinear  model  with  the 
experimental  velocity  records  by  determining  the  parameters 
p  and  k  in  equation  (4.20)  of  the  quasilinear  constitutive 
model.   Numerical  solutions  of  seven  cases  at  prestrains 
equal  to  0.0009,  0.0070,  0.0180,  0.0290,  0.0400,  0.0434  and 
0.0488  were  obtained  first.   The  cost  criterion  in  these 
parameter  identification  studies  was  set  by  comparing  the 
wave  speed  curve  fi  tti  ng  parameters .  v   and  n  of  the  numerical 
solutions  to  the  corresponding  experimental  values  in  each 
case.   The  solution  by  the  quasilinear  model  is  also  not  a 
simple  wave  solution  in  general.   The  maximum  variation  of 
the  propagation  speed  of  a  given  level  of  particle  velocity 
between  values  computed  from  stations  1-2  and  from  stations 
3-4  can,  however,  be  controlled  to  be  under  15%.   The  values 
of  V   and  n  from  the  propagation  speed  between  stations  1-2 
were  used  in  the  parameter  identification  studies.   The  para- 
meter values  resulting  from  these  studies  are  given  in  Table  4.1 
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Figures  3.7  through  3.13  show  comparisons  between  particle 
velocity  records  computed  with  these  resulting  parameter 
values  and  the  experimental  records  at  stations  2,  3,  and  4 
for  the  seven  cases.   Some  remarks  will  be  given  here. 

(a)  The  value  of  n  indicates  the  wave  shape.   This 
factor  is  the  difficult  one  to  match.   It  can  be 
seen  from  the  results  in  Table  4.1  that  the  value 
of  n  from  the  quasilinear  solution  at  e   =  0.0290 

0 

is  much  improved  over  that  of  the  semi  linear  solu- 
tion. 

(b)  The  obtained  parameters  p  and  k  depend  on  the  pre- 
condition parameter  ^  rather  than  depending  on  the 
prestrain  value  alone.  Curves  of  p  and  k  versus  ^ 
are  plotted  in  Figure  4.6 

(c)  It  was  later  found  that  the  choice  of  p  and  k  to 
match  the  particle  velocity  records  did  not  have 
a  unique  result.   Also  an  unreasonable  incremental 
stress-strain  curve  was  found  in  the  results  of 
the  case  e   =  0.0009,  as  will  be  seen  in  Section 
4.5. 

Nevertheless,  at  that  time,  three  other  cases  were 

computed  by  using  values  of  p  and  k  predicted  from  Figure 

4.6.   One  case  is  e   =  0.0121  of  specimen  No.  4  and  the 

other  two  casesare  e   =  0.0100  and  e   =  0.0461  of  specimen 

0  0  '^ 

No.  2.   Cases  of  specimen  No.  2  were  picked  since  the  values 
of  m,  k-j  ,  and  k^  in  Table  3.3  are  close  for  the  two  specimens 
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Results  are  shown  in  Figures  4.7  through  4.9.   It  appears 
that  these  results  are  as  good  as  those  of  the  seven  cases 
previously  shown  in  Figures  3.7  through  3.13. 
{ i  i )   Results  for  aluminum 

Experimental  results  of  specimen  No.  4  were  used  for 
comparison  with  the  numerical  computation.  Eight  cases  at 
prestrain  values  of  0.0015,  0.0060,  0.0100,  0.0140,  0.0200, 
0.0280,  0.0342  and  0.0460  were  examined.  From  the  previous 
experience  in  dealing  with  the  copper  data,  the  quasi  linear 
model  was  tried  first.   It  turned  out  that  for  most  cases 

the  solution  looked  better  if  the  value  of  p  was  on  the  or- 

- 1  2    2 
der  of  10    cm  /dyne  or  pE  -  1.0  (Double  precision  is 

needed  in  the  computation).   It  was  checked  that  these  solu- 
tions were  identical  to  solutions  of  the  semilinear  model. 
Therefore  the  semilinear  model  was  used  to  analyze  the  ex- 
perimental data  of  the  aluminum  specimen. 

Results  of  the  eight  cases  are  shown  in  Figure  3.14 
through  3.21.   It  can  be  seen  from  Figure  3.14  through  3.16 
that  no  reasonable  solutions  were  obtained  for  these  cases. 
For  the  rest  of  the  cases,  the  parameter  k  was  determined 
by  matching  the  velocity  solution  of  the  station  No.  4  with 
the  experimental  data.   From  this  point  of  view,  the  solu- 
tions of  the  semilinear  model  are  only  qualitatively  good. 
The  dependence  of  the  parameter  k  on  the  precondition  factor 
5  is  shown  in  Figure  4.10.   Some  discussion  concerning  the 
incremental  s tress- stra i n  relation  of  the  numerical  solutions 
will  be  given  in  Section  4.5. 
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Table  4.1   Related  Parameters  in  Computer  Solution 
of  Rate-Dependent  Theories. 


Copper 

No.  4 

0.000^: 


0.0070 
0.0180 
0.0290 

0.0400 
0.0434 
0.0488 

0.0121 

Copper 
No.  2 

0.0100 

0.0461 

A 1  u  111  i  n  u  ni 
No.  4 

0.0015   ■ 

0.0060   ■ 

0.0100   ■ 

0.0140 

0.0200 

0.0280 

0.0342 

0.0460 


EXPERIMENTAL 


cm/sec 


0.54   1 .53  0.707 


0.07 
0.39 
0.76 


3.85  0.624 
5.47  0.914 
7.01  0.946 


1 .27  9.43  1 .041 

2.17  14.27  1 . 358 

3.56  18. 71  1.157 

-0. 35  3. 18  0.683 

-0.10  3.35  0. 706 

1 .50  9.81  0.976 


THEORETICAL 


PE 


24.0 
9.6 

16.8 

12.0 


4  - 1 

1 0      sec        cm/sec 


20.9 

18.0 
8.9 


•0.34  0.0 

■0.93  0.0 

■0.18  0.0 

0.16  0.0 

1.67  0.0 

0.76  0.0 

0.50  0.0 

1.68  0.0 


3.5 
6.0 

3.0 

2.0 


10.2 
0.0 

1 
5 

6 

0 

9.6 

0 

5 

6.6 

0 

1 

4.8 
3.0 

0 
1 

1 
0 

3.4 

3.3 
0.6 


3.0 
2.5 
2.1 
1  .8 
1  .5 
1  .0 
1  .2 
0.8 


1  .55 
0.88 


0.68 
0.39 


3.85         0.62 
5.44         0.77 


7.07 
7.66 

9.5 

14.5 

19.6 
22.2 

2.7 


2.83 
10.59 


0.77 
0.56 

0.87 

0.89 

0.92 
0.89 

0.59 


0.57 
0.75 
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Figure  4.6   Relations  Between  Parameters  in  the  Quasilinear 
Model  and  Precondition  Parameter. 
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Figure  4.10   Relations  Between  Parameter  in  Semilinear 
Model  and  Precondition  Parameter. 
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4 . 5   Comparative  Study  of  the  Stra i n-Rate- I ndependent 
and  Strain-Rate-Dependent  Theories 


The  results  for  copper  will  be  discussed  first.   The 
incremental  stress- strai n  curves  obtained  from  the  numerical 
solutions  of  rate-dependent  models  will  be  compared  with  the 
curves  obtained  from  rate- i ndependent  theory.   The  results 
of  the  experimental  data  suggest  strongly  that  a  simple  wave 
solution  exists  in  the  loading  domain  considered.   Also,  the 
curve  fitting  for  the  wave  speed  function  is  a  \/ery    good  ap- 
proximation as  shown  in  Figure  3.24.   The  incremental  stress- 
strain  curvesshown  in  Figure  4.1  are  thus  '^ery    reasonable. 
Examples  of  cases  at  e   =  0.0009,  0.0070,  0.0290  and  0.0488 
are  given  in  Figures  4.11  through  4.14  for  comparison  to 
demonstrate  some  properties  of  the  rate-dependent  models. 

In  Figure  4.11  the  curve  marked  RI  is  the  one  given 
by  the  rate- independent  model,  which  is  believed  reasonable. 
At  the  higher  incremental  strains  it  becomes  parallel  to  the 
quasi  static  curve  indicated  by  the  heavy  solid  line  marked 
QC.   The  curve  marked  QM-Al  is  the  curve  at  station  No.  1 
(the  boundary  for  the  solution),  according  to  the  quasilin- 
ear  model  with  parameter  values  previously  selected  by  para- 
meter identification  to  match  the  velocity  profiles.   It  is 
clear  that  although  this  gives  a  good  match  to  the  velocity 
profiles  (see  Figure  3.7),  the  resulting  incremental  stress- 
strain  curve  is  unreasonable.   It  even  falls  below  the  quasi' 
static  curve. 
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Curve  QM-B]  of  Figure  4.11  resulted  at  station  No.  1 
from  another  solution  with  parameter  choice  B  (pE  =  0.96x10^ 
and  k  =  60,000  sec"  )  .   This  is  a  much  more  reasonable  curve, 
which  lies  close  to  the  rate-independent  curve,  although  its 
slope  at  the  higher  incremental  strains  is  less  than  that  of 
the  rate-independent  curve.   Particle  velocity  predictions 
with  parameter  choice  B  are  shown  in  Figure  4.15.   They  ap- 
pear little  different  from  those  shown  in  Figure  3.7  (ob- 
tained with  parameter  choice  A).   This  shows  that  the  para- 
meter identification  procedure  does  not  produce  unique  re- 
sults for  p  and  k ,  as  was  noted  in  remark  (3)  of  Section 
4.4.5  (i).   The  curves  marked  IC-A  and  IC-B  are  the  instan- 
taneous response  curves  of  the  quasilinear  theory  for  the 
two  parameter  choices. 

This  example,  led  to  the  suggestion  that  the  restric- 
tion of  inquality  (4.38)  on  the  instantaneous  curve  should  be 
strengthened  as  follows  for  use  in  the  parameter  choices: 


dh(c)  ^    E 
de    -  1  +  pEe 


<  E 


(4.64 


where  h(e)  is  the  function  in  equation  (4.6).   Inequality 
(4.64)  is  not  a  necessary  restriction;  it  may  be  violated 
and  solution  may  still  be  a  good  approximation.   From  the 
incremental  stress-strain  curves  of  Figures  4.11  through  4.14 
it  appears  that  this  inequality  (4.64)  is  more  important  for 
cases  with  small  prestrain  or  smaller  values  of  precondition 
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parameter  E,    than  for  cases  with  larger  values  of  E, .       Figure 
4.12  for  the  original  solution  for  the  case  Bq  =  0.0070 
shows  that  a  tendency  similar  to  that  of  the  curve  QM-Al  of 
Figure  4.11  exists  as  the  inequality  (4.64)  is  violated. 
Figure  4.13  gives  a  comparison  between  a  semilinear  solution 
(SM-A)  and  quasi  linear  one  (QM-B).   This  shows  that  the  in- 
cremental stress-strain  curves  of  the  quasilinear  model  are 
more  like  the  ra te- i ndependent  curvesthan  are  those  of  the 
semilinear  model.   Curves  QM-Al  and  QM-A2  in  Figure  4.14 
show  that  a  simple  wave  solution  can  be  approached  as  the 
rate  sensitivity  parameter  k  becomes  small.   The  instantane- 
ous response  curve  IC-A  did  violate  the  inequality  (4.64), 
but  the  result  is  not  so  serious  as  in  the  case  of  Figure 
4.12.   Curves  QM-Bl  and  QM-B2  in  Figure  4.14  show  that  a 
similar  solution  with  a  better  result  in  the  o-e    curve  of 
station  No.  1  is  obtained  by  using  another  pair  of  p,  k 
values.   In  this  case  the  value  of  k  is  larger  and  the  sim- 
ple wave  solution  is  no  longer  approached.   The  particle 
velocity  solution  is  given  in  Figure  4.16. 

From  the  above  analysis,  it  is  concluded  that  the 
relations  between  p,  k  and  E,    shown  in  Figure  4.6  have  to  be 
revised  in  the  region  of  small  E,    and  that  there  exist  many 
pairs  of  the  parameters  p  and  k  which  can  yield  a  good  ap- 
proximate solution.   Values  of  p  and  k  depend  on  precondi- 
tion parameter  C,  but  the  result  is  not  unique.   The  choice 
of  the  form  (1/E  +  pe)  for  (pio,c)    seems  to  be  adequate  for 
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the  small  incremental  waves  in  copper  rods  considered. 

The  inequality  (4.64)  appears  to  be  an  effective 
constraint  to  ensure  that  the  parameters  chosen  to  match 
the  particle  velocity  records  for  the  copper  rods  also  give 
a  reasonable  incremental  stress-strain  curve.   It  also 
proved  to  be  an  expedient  practical  guide  to  a  simplified 
procedure  for  choosing  the  parameters  p  and  k  as  follows. 
First  a  value  of  p  was  chosen    as     large  as  possible  in 
the  expected  c  range,  which  satisfied  the  inequality.   Then 
it  was  found  that  the  parameter  k  could  be  chosen  by  the 
parameter  identification  procedure  to  match  the  velocity 
records.   This  parameter  identification  procedure  was  much 
simplerthan  the  previous  one    of    exploring  combinations 
of  the  two  parameters  p  and  k. 

Results  for  the  aluminum  specimen  will  be  discussed 
next.   In  general  the  solution  by  the  rate-dependent  model 
is  not  satisfactory  quantitatively.   Each  incremental 
stress-strain  curve  in  Figure  4.2  is  almost  bilinear. 
(These  curves  were  obtained  by  using  the  rate-independent 
theory.)   It  is  understandable  from  the  above  analysis  and 
inequality  (4.64)  that  a  semilinear  model  with  instantaneous 
curve  a  =  Ee  is  more  adequate  than  a  quasilinear  model  with 
the  instantaneous  curve  o    =    -   ln(l  +  pEe).   By  checking 
through  Figures  3.14  to  3.21  it  is  seen  that  in  general  the 
solution  by  a  semilinear  model  is  qualitatively  satisfactory 
for  most  cases  except  the  cases  at  e^  ^  0.0015,  0.0060  and 


131 


0.0100.   One  important  fact  supporting  the  solution  of  the 
semilinear  model  is  that  in  the  experimental  records  the 
particle  velocity  plateau  approached  decayed  along  the  prop- 
agation  distance  x.   This  dynamic  relaxation  phenomenon 
cannot  be  explained  by  the  simple  wave  solution  of  rate- 
independent  theory.   Most  incremental  stress-strain  curves 
of  Figure  4.2  were  obtained  by  assuming  that  the  particle 
velocity  continuously  increased  to  40  cm/sec,  which  was  not 
true  in  the  experiments.   As  shown  in  Figure  3.23,  many  wave 
speed  functions  ceased  at  a  certain  velocity  level.   An  exam- 
ple of  an  incremental  stress-strain  curve  obtained  from  the 
numerical  solution  with  the  semilinear  model  is  given  in 
Figure  4.17  for  comparison  with  the  rate- i ndependent  one. 
The  stress  relaxation  shown  in  the  semilinear  solutions  is 
qualitatively  not  unreasonable  in  view  of  the  continued 
straining  at  the  particle  velocity  plateaus. 

About  the  results  of  aluminum,  more  study  is  needed. 
From  the  results  of  the  experiment  it  appears  that  the  be- 
havior of  aluminum  is  very    complex.   We  can  conclude  only 
that  the  small  incremental  wave  can  be  described  qualita- 
tively by  using  a  semilinear  model  of  the  Malvern  type  for 

cases  at  e„  >  0.0100.   The  rate-sensitivity  constant  k  is 
0  -^ 

also  found  to  be  dependent  on  the  precondition  parameter  F, . 


CHAPTER  5 
SUMMARY  AND  CONCLUSIONS 

A  series  of  incremental  wave  tests  were  performed  on 
annealed  aluminum  1100-F  and  commercially  purecopper  rods 
under  quasistatic  preloading.   An  electromagnetic  transducer 
was  used  to  obtain  measurements  of  particle  velocity  at  four 

stations  along  the  rod.   The  strain  rates  of  the  quasistatic 

-6      -4    -1 
preloading  were  in  the  range  of  10    to  10    sec   ,  while 

strain  rates  of  the  order  of  1.0  sec   were  achieved  in  the 

i  ncremental  waves . 

Test  results  indicated  that  except  near  the  plateau 

in  the  records  for  aluminum  any  given  level  of  incremental 
particle  velocity  propagates  along  the  rod  with  a  constant 
speed,  not  affected  by  the  strain  rate  within  the  small 
range  of  the  incremental-wave  strain  rates  encountered. 
However,  the  wave-speed  function  of  particle  velocity  dif- 
fered noticeably  for  the  incremental  waves  generated  at  dif- 
ferent prestrain  levels.   An  exponential  function  with  two 
parameters  was  found  to  be  able  to  fit  the  various  wave 
speed  functions  of  both  materials.   The  results  then  in- 
dicated that  one  of  the  parameters,  the  major  one,  v,  de- 
pends not  only  on  the  level  of  prestrain  but  also  on  the 
value  of  prestrain  rate.   This  discovery  shows  another  phy- 
sical property  of  metal  under  incremental  loading,  in 
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addition  to  the  well-known  phenomenon  of  the  leading  edge 
propagating  at  the  elastic  bar  wave  speed. 

This  dependence  on  the  prestrain  rate  may  be  related 
to  a  phenomenon  observed  by  Klepaczko  et  al .  (1974),  who 
found  the  response  to  a  sudden  rate  change  in  a  torsional 
test  to  be  dependent  on  the  preceding  strain  rate  history. 
They  did  not  examine  wave  propagation  but  measured  average 
stress-strain  response  along  a  short  specimen. 

Some  preliminary  efforts  were  made  to  describe  quan- 
titatively the  dependence  of  the  incremental  waves  on  both 
prestrain  and  prestrain  rate.   A  new  parameter  composed  of 
these  two  variables  was  introduced,  which  was  called  the 
precondition  parameter  ^  in  the  text.   The  possibility  of 
interpreting  the  pres tra i n- ra te  dependence  from  the  point  of 
view  of  the  strain-rate  difference  between  the  incremental 
wave  and  quasistatic  preloading  should  not  be  excluded,  al- 
though not  enough  evidence  exists  to  show  this. 

Both  strain-rate-dependent  theory  and  strain-rate- 
independent  theories  were  compared  with  the  results  of  the 
experiment.   Based  on  the  simple  wave  solution  of  rate- 
independent  theory,  incremental  stress-strain  curves  were 
obtained  for  both  cases  of  aluminum  and  copper.   Quasilinear 
and  semilinear  rate-type  constitutive  models  with  linearized 
terms  of  ())(a,e)  and  \])io,e)    were  proposed.   Numerical  com- 
puter solutions  were  obtained  for  incremental  waves  at  vari- 
ous prestrain  levels  of  both  aluminum  and  copper  specimens. 
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The  velocity  records  obtained  from  the  first  gage  station 
(15  diameters  from  the  impact  transmitting  collar)  were  used 
as  input  boundary  conditions. 

An  i ntegro-di f ferential  approach  was  introduced  for 
obtaining  numerical  computer  solutions  with  the  quasilinear 
model.   It  turned  out  that  this  technique  is  a  very   effi- 
cient one  in  solving  the  system  of  partial  differential 
equations  of  nonlinear  wave  propagation. 

Numerical  solutions  show  that  the  proposed  quasilin- 
ear constitutive  equation  is  an  adequate  one  to  describe  the 
copper  material.   In  the  copper  cases,  it  was  found  that 
matching  the  particle  velocity  records  does  not  guarantee 
correctness  of  the  solution.   An  additional  check  on  the 
solution  was  made  by  comparing  the  incremental  stress-strain 
curve  with  the  stress-strain  curve  obtained  from  the  rate- 
independent  theory.   When  the  two  agreed  approximately  a 
good  match  to  the  particle  velocity  records  was  always  ob- 
tained.  Therefore  a  stronger  inequality  condition  was  pro- 
posed [inequality  (4.64)],  which  proved  to  be  a  sufficient 
(though  not  always  necessary)  condition  for  determining  the 
correct  range  for  parameter  p  in  the  instantaneous  plastic 
term  (j)(a,e).   By  adjusting  the  two  parameters  in  the  quasi- 
linear  constitutive  equation,  a  simple  wave  solution  with 
its  different  dynamic  stress-strain  curve  for  each  prestrain 
level  could  be  well  approximated  by  the  quasilinear  model. 
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The  experimental  results  for  the  aluminum  specimen 
show  the  phenomenon  of  dynamic  relaxation  which  could  not 
be  described  from  the  simple  wave  solution  of  the  strain- 
rate-  i  ndependent  theory.  Computer  solution  with  the  semi- 
linear  model  was  able  to  describe  it  qualitatively.  However 
more  efforts  should  be  made  in  the  future  to  describe  the 
experimental  data  of  incremental  waves  generated  at  small 
pres tra  i  n  1 evel s . 

Further  study  in  this  area  should  include  the  fol- 
lowing:  (1)  A  series  of  tests  with  increasing  amplitude 
of  impact  for  generating  the  incremental  waves  would  be 
useful,  since  larger  amplitude  usually  results  in  larger 
values  of  strain  rate.   (2)  The  quasistatic  preloading 
should  be  controlled  at  a  given  constant  prestrain  rate  for 
a  test  of  one  speciman.   Incremental  wave  tests  should  be 
performed  on  several  specimens  at  various  constant  prestrain 

rate  values  differing  from  each  other  by  a  magnitude  of  the 

-5   - 1 
orderof  10  sec   .   (3)  A  longer  copper  rod  is  preferred  so 

that  reflection  from  the  end  will  not  reach  the  gage  sta- 
tions before  particle  velocity  reaches   a   constant  level. 
The  dynamic  relaxation  phenomenon  can  then  be  checked  in 
copper. 


APPENDIX 


COMPUTER  PROGRAM 


A-1   Velocity  Data  Processing  Program 
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